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VECTOR CALCULUS 




INTRODUCTION 

In course of an attempt to apply direct vector methods to 
eeitaia problems of Electricity and Hydrodviiamics, it was 
felt that, at least as a matter of consistency, the foundations 
of Veetoi Analysis oug'ht to be placed on a basis independent 
of any reference to cartesian coordinates and the main theorems 
of that Analysis established directly from first principles. 
The result of my work in this connection is embodied in the 
present paper and an attempt is made here to develop the 
Differenhal and Integral Calculus of Vectors from a point of 
view which is believed to be new. 

la order to realise the special features of my presentation 
o the subject, it will be convenient to recall briefly the usual 
method of treatment. In any vector problem we are <riven 
certain relations among a number of vectors aud we "have 
to deduce some other relations which these same vectors 
satisty. Now what we do in the usual method is to resolve 
each vector into three arbitrary components and thus rob it 
hrst entirely of its vectorial character. The various charaoteris- 
le vector operators like the gradient and curl are also subjected 
to the same process of dissection. We then work the whole 
pro em out with our familiar scalar calculus, and when the 
necessaiy analysis has been completed, we collect our components 
an re the result in vector language. It is of course quite 
e u so ai as it goes, the final vector expression of the result 
giving not only a succinct look to our formula; but also a 
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suggestiveiiess of interpretation which they had been lacking 
in their bulky cartesian forms. But surely, strictly speaking, 
we should not call it Vector Analysis at all, but only Cartesian 
Analysis in vector language. In Vector Analysis proper we 
have, or ought to have, the vector physical magnitudes which 
our vectors represent, direct before our minds, and this charac- 
teristic advantage of being in direct close touch with the only 
relevant elements of our problem is sacrificed straight away, 
if we throw over our vectors at the very outset and work with 
cartesian components. We sacrifice in fact the very soul of 
Vector Analysis and what remains amounts practically to a 
system of abridged notation for certain complicated formula 
and operators of cartesian calculus which happen to recur every 
now and then in physical applications. 

The one great fact in favour of this plan is that it affords 
us greater facility for working purposes, this facility no doubt 
arising solely from our previous exclusive familiarity with Carte- 
sian Analysis. But however useful it might be in this direction, 
and generally in making the existing body of Cartesian Analysis 
available for vector purposes, the process, I venture to think, 
is at best transitional, and the importance of the subject and 
the importance of our thinking of vector physical magnitudes 
direct as vectors, alike seem to demand that the whole of this 
branch of Analysis should be placed on an independent basis. 

But there is a peculiar difficulty at the very outset. Histori- 
cally, most of the characteristic concepts of Vector Analysis, 
like the divergence and curl, had been arrived at by the physicist 
and the mathematician in course of their work with the Cartesian 
calculus and had even become quite familiar before the 
possibility of Vector Analysis as a distinct branch of mathematics 
by itself was explicitly recognised. The vector analyst at 
first then starts from these old concepts which happen also 
to be the most fundamental, but it is his object right from the 
beginning to exhibit them no longer in their cartesian forms, 
but in terms of the characteristic physical or geometrical 
attributes which they stand for. Very often now a question 
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of selection arises from among the number of ways in which 
the saime concept may be defined, different definitions being 
framed according to the different j)oints of view from which 
the subject is intended to be developed. The physicist — who, 
by the way, makes the greatest practical use of Vector Analysis 
and whose sole interest also in the subject is determined by 
the service it renders him in his work — aims, first of all, at 
his definition representing most directly a familiar physical 
fact or idea; but, at the same time, and very naturally too, he 
holds the possibility of the definition yielding quite easily 
his useful working formulie, of equally vital importance. But, 
unfoilunately enough, these two distinct aims of the physicist 
are irreconcileable with each other, the most natural definition 
from the physical point of view leads to the useful transfor- 


mation formulse of Physics only with the greatest difficulty, 
and the definition that yields these formulae with any facility 
IS generally hopelessly artificial from the physical point of 
view.* It is this irreconcileability of the two distinct purposes 
of the physicist which, I venture to suppose, is directly respon- 
sible for the persistence of carte.sian calculus in Vector Analysis. 
For what is done is that definitions are first framed with a 
view to direct summing up of the simplest appropriate 
physical ideas, but then the necessity almost inevitably arises of 
seeking cartesian expressions for working purposes, for making 
Vector Analysis a serviceable and at the same time an easily 
manageable tool in the hands of the physicist. 

I may just illustrate my point by recalling how the usual 
definitions of the two most characteristic concepts of Vector 

. Koference may' be made here to a paper by in the 

^l eHn of the Amcr^an Math. Hoc, vol, 16, on UniBcation of Vectorial 
Notatmae, where he eritieieeB the artificiality in the definitions of divergence 
and curl by an Italian mathematician, Burali Forti, which were chosen 
Bolely with a view to their adaptability for establishing the working formula 
of Vector Analysis with ease. Thus Burali Forti’s definition of divergence is 
div V = a.[grad (y) + ourl (axv)], where a is any constant unit vector 
This has certainly no direct connection with any intrinsic property of the 
divoi-genc®, phyaioal or otherwise. t' k> y oi me 
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Analysis have been adopted from the simplest physical ideas 
which immediately identify them. Thus the idea of divergence 
is taken directly from Hydrodynamics, and keeping before our 
minds the picture of fluid leaving (or entering) a small closed 
space, we define the divergence of a vector function at a point 
as the limit of the ratio, if one exists, of the suifaee-integi'al * 
of the function over a small closed space surrounding the point 
to the volume enclosed by the surface, a unique limit being 
supposed to be reached by the closed surface shrinking up to a 
point in any manner. 

Again, it is found that some vector fields can be specified 
completely by the gradient of a scalar function, so that the line 
integral t of the vector function along any closed curve in 
(simply connected) space would vanish. Thus the work done is * 
nil along any closed path in a conservative field of force. But 
in case the vector function cannot be so specified, an expression 
of this negative quality of the function at a point is naturally 
sought in its now non-evanescent line integral along a small 
closed (plane) path surrounding the point. The ratio of this 
line integral to the area enclosed by our path generally approa- 
ches a limit as the path shrinks up to a point, independently of 
its original form and of the manner of its shrinking, but dep- 
ending on the orientation of its plane. The limit moreover has 
usually a maximum value, subject to the variation of this 
orientation, and a vector of magnitude equal to this maximum 
value and drawn perpendicular to that aspect of the plane which 
gives us the maximum value is called the curl of the original 
vector function. 

Now these definitions, embodying, as they do, the most 
essential physical attributes of divergence and curl, must be 
regarded as perhaps the most appropriate ones that could be 
given from the physicist’s point of view. But then comes the 

* By the surface integral of a vector function, we always mean the surface 
integral of its normal component. 

t By the line integral of a vector function, we always mean the line inte- 
gral of its tangential component. 
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practical problem of deducing from these dehnitious the working 
rules of manipulation of these operators. The direct deduction 
being extremely difficult,* the already acquired facility in 
working cartesian calculus is naturally utilised for the purpose, 
and thus is reached the present position of Vector Analysis 
which I have already described. 

The only way out of the dilemma would seem to be found 
by ignoring altogether both of these two specific interests of the 
physicist and looking straight, without any bias, to the require- 
ments of Vector Analysis as a branch of Pure Mathematics bv 
itself. And paradoxical though it may sound, this course 
perhaps would ultimately best serve the physicist’s ends also. 
At any i-ate, no free development of any science is certainly 

possible, so long as we require it at every step to serve some 
narrow specific end. 


We ask ourselves then, what should be the most natural 

starting point of the DifEerential Calculus of Vectors ? All our 

old familiar ideas of differential calculus suggest at once that, 

whatever the ultimately fundamental concepts might be, we 

should begin by an examination of the relation between 

the differential of the vector function (of the position of a point 

P m space) corresponding to a small displacement of the 

point P and this displacement. This very straightforward line 

of enquiry I propose to conduct here, and it will be seen how 

in yery natural sense we can look upon the divergence and 

curl M really the fundamental concepts of the Differential 

Calculus of vectors, and how this new point of view materiallv 
Simplifies our analysis. 


The first three sections are preliminary. In the first two I 

summarise the definitions of continuous functions and of Inte- 
gra s and briefly touch upon just those properties which I 
require in course of my work. The third is devoted to the 
Gradient of a scalar function. T he real thesis of the paper I 

encountered by Mr. E. Cunning- 
of the'ZS^wLT Proceeding. 
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begin in the fourth section where 1 consider the Linear Vector 
Function only with a view to developing what I have called the 
scalar and vector constants of the linear function, and although 
there is nothing very special about these ideas themselves, they 
will be found to lead very naturally to the concepts of Diver- 
gence and Curl and have been made here the foundation on 
which my Diffei*ential Calculus is built. The fifth section is 
devoted to that Differential Calculus and in the sixth I consider 
a few Integration Theorems and the divergence and curl of an 
integral with a view to showing with what ease these operations 
can be performed from ray point of view. 

Notation. 

With regard to notation I use Gibbs’ here, although some 
of its features are obviously meant to suggest easy ways of 
passing from Cartesian formulae to vector, and vice verm, with 
which of course I am not at all concerned. 

For convenience of reference I reproduce the notation for 
the multiplication of vectors. 

If A, B are any two vectors, 

the scalar product of A, B is A,B= | A | | B | cos Q, and the 
vector product is A x B which is a vector of magnitude | A | 

I B I sin and in direction perpendicular to both A and B ; 

I A I , I B I denoting the tensors of A and B and B the angle 
between them. 

Again, if A, B, C are any three vectors, the notation 
[ ABC ] is used for any one of the three equal products 

A.B X C = B.C X A = C.A X B = the volume of the paiallelo- 
piped which has A, B, C, for conterminous edges. 

The following useful formula will occur very often : 

A X (B X C) = (A.C)B-(A.B)C. 



I. 


Continuity : Differentiation of a Vector 
Function of a Scalar Variable. 


1 . The functions we deal with will be mostly continuous. The 

position of a point P in space being specified as usual by the 
vector »’( = OP) drawn froir. a fixed origin O, the function f{r) 
IS said to be continuous at P, if corresponding to every arbitra- 
rily chosen positive number 8, a positive number r, (dependent 
on 8) can be found such that /(r-i-t) _/(,•) , <8, c being any 
vector satisfying the inequality i € i < The notation iVi 

denotes the absolute value of the scalar if V is a scalar, and the 
tensor of V if V is a vector. 


If we construct the vector diagram as well, that is, if by 
taking another fixed point O' we dmw the vector O'P' repre- 
senting the value of the vector function corresponding to every 
point P in the region in which the function is defined, then Q 

I ^ of P and Q' the correspond- 

ing point in the vector diagram, our definition of continuity 

implies that any positive number 8 being fii-st assigned, a 

positive number , can be found such that so long as the tensor 

of the vector PQ, is less than -q, the tensor of P'Q' will be less 

than 8. It implies in other words that a sphere (of radius q) 

can be described with centre P such that points Q' in the vector 

diagram corresponding to all points Q within (not c«) this 

sphere will he within a sphere of any arbitrarily small radius 
S described with centre P'. 

We prove now that in the same case the angle P'O'Q' 
at IS, the change in direction suffered by the vector fune- 

mn Mn also be made arbitrarily small. For, in the triangle 
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A A 

P'O'Q' _ sin P'Q'O' ^ 1 

P'Q' 

Hence, sin PU'Q' F Q 

^ ^ O'P' 

But P Q can be made arbitrarily small, and O'P' is 
supposed to be finite. Hence sin P'O^Q' and therefore also the 
angle P O'Q,' can be made arbitrarily small. It follows that our 
continuous vector functions are continuous in direction as well. 

The function f (?•) is said to have a limit at P, if Q 
being any point in the neighbourhood of P we have the same 
limiting value of the function no matter in what maimer Q 
approaches P continuously. 

If the function is continuous at P, the limit exists at P and 
is e<^ual to the value of the function at P, and conversely. 

If the limit does not exist at P, then either of two things 
may happen : [t) there may be different limiting values for 
different approaches to P ; or (//) there may be no definite 
limiting value for any approach or some approaches. In either 
case the function is discontinuous at P. 

A third kind of discontinuity arises when the limit exists 
at P, but this limit is not equal to the value of the function at P. 

But, as has been remarked already, we shall concern our- 
selves practically with continuous functions alone, and an 
examination of the sort of peculiarities we have just noticed, of 
what has been described as the Pathology of Functions would 
be out of place here. The only discontinuity we shall come 
across is the infinite discontinuity which arises when | 

tends to infinity at P. 

3. Turning our attention then to continuous functions 
alone, we note that the sum and the scalar and vector products 
of two continuous vector functions are continuous also. The 
ease of sum is almost self evident, and we prove now that if V,, 
Vg are two continuous vector functions, the scalar product 
V,-Va is continuous. 
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Let V\, denote the values of the functions at a point 
in the neighbourhood of the point We have only to show 
that for any positive number 8 assigned in advance, a positive 
number ly can be found such that 

1 1 < 8 , 

for all vectors c satisfying ] c | <?y. 

Now + (V\-V,)]-[V, 

+ (W.-VJ]-V, ■V,=V, -(V.'-VJ 

which is not gi'eater than j V, | | Vj'— Vjj | + | V 2 | | Vj-V, | 

+ I Y/— Vj 1 I Vg'— Vj I , since the magnitude of the scalar 

product of two vectors is not greater than the product of their 
tensors. 

HencOi since the absolute magnitude of the sum of any 
number of quantities is not greater than the sum of their abso- 
lute magnitudes, we have 

lV/-V/-V,-Vg|>lVJ IVg'-Vgl 

+ 1 VJ I V.'-V, I + I V,'-V, 1 I V,'-V, I 
But since V,, V, are continuous, 

I ^\~y \ 1 <any arbitrary 8,, provided only | € | < the 

corresponding ,, and | V,'-V, | < 8, -q^ 

Of the two numbers i;, r;,, let ; 

then provided 1 t \ \ V/-V, \ <8. and 1 V,'-V, | <8„ 

and therefore 1 V,'- V,'-V, • V, | < 1 V, j 8, + 1 V, | 8, +8,8,. 

Again, since | V, | , | V, | are supposed to be finite, given any 
positive number 8, we can always find 8^ and 8g such that 8> 

18.+ I Vg 18,+8,8g. 

Choosing such values now of Sj and 8,, we have 

I V\*V\— 'Vi Vg I <8, whenever | e | < 1 /^, 
which proves our theorem. 

Similarly we prove that Vj x Vg is also continuous. 
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4. If we consider in particular the continuous vector func- 

tion of a scalar variable, we can easily adapt the argument of 
the usual scalar calculus and prove the theorems associated with 
continuity in that calculus. If v—f{t) be the function consi- 
dered, t being the scalar variable, we can prove in particular 
that if Tx—f (/^i) and (^ 2 ) P ^”7 number 

such that I t*! I <p< I ^*2 I > there is a value of t lying 

between and for which |,/(/) | =p. In other words, as t 
varies continuously from to the tensor of v assumes at least 
once every value lying between and 

It can further be proved that if F (^) is any continuous func- 
tion, scalar or vector, of v where v itself is a continuous function 
of a scalar variable t, then F is a continuous function of t * In 
case F is a scalar function, it follows that if Fj, are 
the values of F respectively for ^ = and t — then as t 
varies continuously from to ^ assumes at least once 
every value lying between F^ and Fg ; and when F is a vector 
function, it is the tensor of F that assumes, as t varies conti- 
nuously from ix to least once every value lying between 

the tensors of F corresponding to ^ = and 

5. If for the continuous vector function a unique 

limit exists of ^ as t' approaches t from either side, 

t'-t 

from values less than ^ to ^ and from values greater than 
t to 1), then this limit is called the difFerential coeflScient of 

* Proof. We have to show that if 5 is assigned in advance, 1? can be 
found such that 

IFCO-FO*) 1 

when 1 f'— < 1 AVt being the value of r corresponding to t=t\ 

Now since F(r) is a continuous function of r, an 17, can be found such that 
1 F(r') — F(r) 1 ^5, when 1 / — r 1 Zi/,. 

Again, since r ia a continuous function of t, corresponding to this 9711 & 
positive number 17 can be found such that 1 /— r 1 AVa when 1 1' — 1 1 Av- 

This 77 then is such that when 1 t'—t 1 Av, 

1 r'— r 1 AVi 1 — F(r) 1 ^77, 
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i* with respect to i and is denoted by 



The function in 


the same ease is said to be differentiable at 


A function =/’(0 which is continuous and differentiable 
at all points in a certain region can in general be represented 
by a curve in that region. The terminus of v will trace out 
the curve as t goes on varying continuously^ and the vector 



will be at each point in the direction of the tangent to 


the curve at that point. 

If F(j*) is a continuous vector function of it follows now 
fiom the last article, that the vector diagram of F(i*) eorres- 
l>onding to points lying on any arbitrary but continuous curve 
r=^(f) between any two specified points P and Q is also a 
continuous curve lying between the corresponding points P' 
and Q,' in the vector diagram, 

6. Mean value theorem for — If y is a continuous 

and differentiable function of t for all values of t between any 
two specified numbers f ^ and then and being the 
values of r respectively for ^ = and / = ^ 2 , we have 


where B is some positive proper fraction. 

This IS proved, precisely as in the case of the corresponding 
theorem in scalar calculus, by considering the function 




• 2 




^2 




which is continuous for all values of / between t, and and 

nishes foi t and and of which therefore the differ- 

ential co-efficient will vanish at some point between and 

say at where ^ is a positive proper fraction. This 

proves our theorem. 
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Graphically, if R denotes the vector to any point on the 
chord of the curve v=f(t) joining the points and the 
equation of the chord is 


R=»'i 

for obviously it represents for varying values I a straight line 
parallel to — and gives R=:a'j at ^ = and R = r2 at 

Our = R = R', say, represents then for any value 

of I the difference of the vectors to points on the curve and 
the chord corresponding to that value of t. If these vector 
differences are now drawn from the origin for all values of i 
from to /gj their terminii will give us another curve re- 
presented by R' = 4 >(/), which clearly is a continuous curve 
returning unto itself at the origin for ^=^1 and and the 


vanishing of 


dt 


for some intermediate value of t implies that 


in course of the journey of the terminus of R' from the origin 
and back to it again, there will be a position which will make 
the tensor of R' or r— R stationary. 



II. 

Integrals. 


7. The Fedor Volume Integral . — Given any finite con- 
tinuous volume T, if for any eonverg-ent system of sub- 
divisions * of the region, the vector sum ^ where 

T, denotes any sub-region at any stage of the sub-division 
and F, the value of F at any point within the sub-region 
T., tends to a definite, unicjue limit as the sub-division 
advances, independent of the particular convergent system of 
sub-d.v.sions used and of the particular values of F chosen 
within the sub-regions r.. then this limit is called the volume 
integral of the vector function F through the volume r, and 

T 


is written J Ydr. 

Without gomg into the question of the necessary minimum 
condition for the integrability of F, we may prL wiITt 
much trouble the only theorem we require in this connec- 
lon, VIZ., that if F is continuous at all points within a 
e reyon, it is integrable also through that region — 
e continuity of F ensuring that if F., F/ are the values 

r sub-region r„ the tensor of 

the difference of P. and F.' becomes arbiti-arily small as 

volume sub-region diminishes in 


mav aSo J ‘'^Presentation of the vector volume integral 

pont O, we lay down the vectors F.r. as in the ordinary 
^lygon of vectors. In the limit the polygon becomes I 
continuous curve, ending say in A.' Then the arc O'A' will 

represent / | P | rfr and the chord oa' will represent our 
volume in tegral /Frfr. Since F is supposed to be integrable, 

* Compare Hobaon’e Theory of F^^Tuone of a Real Variable. § 251. 
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the chord O'A' will be unique, but we may have an infinite 
number ol curves like O'A' according to the different orders in 
which we may place the vectors F„t„ in forming the 
polygon. All these curves however will have the same 
length / I F I /'/r and the same terminal point A'. Further, to 
any point P' on any one of these curves there will corres- 
pond a unique point P in the volume t, and conversely, 
so that there is a one-to-one correspondence between the points in 
the volume and the points on any particular curve.* We 
have also fh*' = V'fhf if is the vector O^P', so that the 
tangent at any point P' on the curve is in the direction 
of F at the corresponding point in the volume. It may 
happen that / | F | (It is infinite, but JF^/t at the same time 
exists as a finite vector. Thus the curve may make an 
infinite number of convolutions, but such that the terminal 
point A' is at a finite distance from O'. 

8. T//e surface integrals , — Given any continuous surface 
S in a region where the vector function F is defined, we 

o 

form the scalar product F.m at each point of the surface, 
n denoting the unit vector along the outward normal at 
any point to the surface, and the surface integral, in thQ 
usual sense, of the scalar function F.m over the surface we 
call the surface integral of vector function F over the 
surface and denote it by /F.«^^S. In other words, for any 
convergent system of sub-divisions of the surface S, if Sr 
is a sub-area at any stage of the sub-division and 
the value of Y.n at any point within the sub-area, the unique 
limit to which is assumed to tend as the sub- 

division advances is called the suiface integral of F over 
the surface S. But wdth the advance of the sub-division 
the areas S, approximate to small plane areas on the 
tangent planes at points P, and the vectors ?^rST ultimately 

* There is of course no a priori absurdity in the idea of a one-to- 
one correspondence being established between the points in a given 
volume and the points on a line, for we know from the theory of 
sets of points that the two aggregates have the same '^power.” 
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may be regarded as representing these plane areas both in 
magnitude and direction. A\e may replace tlierefore the 
notation jF.w^/S iby //cr=r//^/S representing the ulti- 

mately plane element ^/S both in direction and ma«’nitnde. 

© 

Forming again the vector product of F and the vector 
at each point and summing np for all j)oints and 
passing to tlie limit in the same wav, we have another 
surface integral {FxnifS or /Fx^/<r. This has been called 
the skew or vector surface integral, /F.^/o- being the direct 
or scalar surface integral. We shall always mean jF.^/o- 

when we speak only of the surface integral of F, referring 
to JF X dtr as the skew surface integral. 

o 

If the surface S is supposed to be continuous and to possess 

moreover a continuous tangent plane at every point, the vector 

n would be a continuous function over the surface, and if F 

IS supposed to be continuous also, both F. n and Fx n will be 

continuous functions and the scalar and vector surface inteo-rals 

of F oyer S will both exist. We shall always mahe”this 
supposition here. 

9. Tke lute tnteffralfi.—G\\e.n any continuous curve, if in 
any convergent system of subdivisions, is the vector chord 
.loimng two consecutive points of division at any stage of sub- 
division in the.system, and F, is the value of F at any point 

P of the curve between the.se two points of division, then the 
unique limit to which 5F..P. is assumed to tend as the .sub- 
division advances is called the line integral of F alon- the 
curve AB The chord p, is obviously equal to the difference 
in the values of r at the two points which it connects, and our 

integral may be denoted by It is further clear that 

with the advance of the subdivision, p, approaches in direction 

to the tangent to the curve at P, and if therefore we denote 

he unit vector along the tangent at any point of the curve by 

t, the integral is the same as the line integral, in the usua^ 
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sense, of the scalar function F.t, and might be denoted by / F.t ds, 
ds being the scalar element of arc. 

We might in the same way define the vector line integral 
/ F X dr, but this wdll rarely occur in the present paper. 

In any ease we shall always suppose that the curve along 
which we integrate is not only continuous, but also possesses 
a continuous tangent, so that t is a continuous vector function 
of the position of a point on the curve. 

The following properties of the line integral follow imme- 
diately from the definition. 



rB CP fH 

[ii) 1 F.dr=\ F. di'H- I F.dr, P being any point on the 

)a ja jp 


curve AB. 

{iii) If 1 is the length of the arc AB and L, U the lower 
and upper limits respectively of F.t for the curve AB (which 
limits are supposed to exist, though not necessarily to be 

attained), then 


Ll^ 


rB 

I F. dr 


Ul. 


(^?^) Further, if M is some number satisfying LrSM^U, 


we have 


TB 

J F.dr = Ml ; 


and in ease F is continuous, so that F.t 


IS 


continuous also, the value M is attained by F.t at some point 


P of the curve, and we have 


rB 

1 F.dr = (P.t)=l. 


Ill 


The Gradient of a Scalar Function. 


10, Let F(r) be a continuous scalar function of the position 
of a point P (OP = r) in a given region. If Q is a point in the 
neighbourhood of P, such that PQ = ah where a is a unit vector in 
direction PQ and h a small positive number, then the value of 

F at Q is P (r f ah). If now the limit [F (r H- ah) — 

^(0] exists as a definite scalar function (difPerent from zero) 
of a and r, this limit would measure the rate of change in the 
value of the function for a displacement of P in the direction 
a* Supposing the limit to exist and denoting it by f (a, r), 

we have F(rH-ah) — F(r)=hf(a,v)-|-hi7, where and h have 
the simultaneous limit zero. 


Now since h appears in the left hand side of this equation 
only in the combination ah, and the first term on the right hand 
side is linear in h, it follows that this term is linear in a also. 


Ihe function f (a, r) then is a scalar function linear in a ; it 
vanishes moreover with a, and therefore it must be of the form 


a*- G (r), where G (r) is a vector function of r, independent 

of a. 


If the limit in question exists now for every direction a 

emerging from P, the i-ate of change of F(r) in any direction 

a is a.G(r), the maximum value of which obviously, for varying 

directions a, is obtained when a is taken in the direction of G 

and the magnitude of the maximum value is equal to the tensor 

of G. ^ The vector G is called the gradient of the scalar function 

P. The gradient of a scalar function then may be genei-ally 

defined as a vector in the direction of the most rapid i-ate of 

increase of the function and equal in magnitude to this most 
rapid rate. 

3 
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il. The same question may be looked at geometrically also. 
We begin by proving* that if F (r) is a scalar function continu- 
ous in a certain region and does not possess any maxima or 
minima in the region, and if is the value of the function at 

any point P of the region, then there passes through P-a surface 
on every point of which F has the value F 

^ 9 

For, since P is neither a point of maximum nor miniraum, 
all the values of F in the neighbourhood of P cannot be greater 
than Fp , nor can all the values be less than F^ , and there would 

be points in the neighbourhood for which F is greater than Fp 

and there would be points also for which F is less than Fp . 

In the neighbourhood of P then, let Q be a point such that 

F >F , and R a point such that F <F„ • Now on account 
Q P R P 

of the continuity of the function, a region can be constructed 

about Q within which the fluctuation of the function is as small 

as we please. Hence there exist other points near Q for which 

also the value of the function is greater than Fp . Similarly 

there exist points near R for which the function is less than Fp • 

Hence the region consists of two distinct regions in every point 
of one of which F>Fp , and in every point of the other 


F<Fp. 

Again, since in passing from any Q to any R along a conti- 
nuous curve, F must on account of its continuity assume all the 
intermediate values, it assumes the value Fp somewhere between 


* This proof is adapted from the eolation an example in Ronth's staticSi 
Vol. II (Ex. 2, $ 124), where from the fact that gravitational potential is 
neither a maximum nor a minimum in free space is deduced that an isolated 
line cannot from part of a level surface. 
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Q and R on that curve. Hence there is a continuous surface 
of separation of the two regions at every point of which F=P^ , 

which proves our theorem. 


If now the surface possesses a tangent plane at P, we take 
a point P' on the normal to the surface at P in its neishbourhood. 
Through this point P' also will pass a surface on every point 
on which F=F ^ and PP' will be normal to both the surfaces 


F=Fpand F^F^, 


Supposing now that the limit 


F^^ - F 

ppr- 


exists as P moves continuously along the normal and approaches 
P, a vector iu the direction of this normal and equal in magni- 
tude to the value of this limit is called the gradient of F(;*) at P. 
[P might be on the normal on either side of the surface, and 
it is assumed that the limit in question exists in either case and 
that these two limits are equal.] 

To see that the gradient so defined gives us the most rapid 

rate of increase of the function both in magnitide and direction, 

we take a point Q in the neighbourhood of P on the surface 

on which P lies. Let ZP'PQ=:d. Then the rate of increase 
of the function in direction PQ. 

■p p 

T Q ^ T F F PP' 

=QP=,0 ~PQ PF=0 ~PP^' 

of which the maximum value obviously is obtained when 6=0, 

This establishes the identity of the definitions of gradient in 
the present article and the last. 

We denote the gradient by VP. If 8P is the change in the 
value of P on account of the shift Br in the position of P, we have 

8F=VP. Br+r, I Sr I 

where and \ Sr | have the similtaneous limit zero. 

12. If the shift Sr is supposed to take place along a definite 
continuous curve r=x(t), then as we have seen (§4, p. 1 0) F 

would be a continuous function of t along that curve, and our 

% 

relation of the last article can be written vF. 

dt di * 


no 
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Further, if sj^ecify any two points K, L on the curve 
x(0 if ^ 1 ) Fg are the values of F at K and L respec- 
tively, we have by the Mean Value Theorem of §6, p.ll. 



where dennotes the value of at some point M on the 

dt ^ 

curve lying between K and L. That is to say, 


F.-P, (VP)., 

In particular, if the curve is a straight line in the direction of 
the (unit) vector a and k is the length of KL, so that 
we may write F(r-|-rt/4)— F(/-)=^fl.(VF)«=//a. vF 
where ^ is a positive proper fraction ; or again, P (r+a)— F(^’) 
=a.vF(';-H-^a), 

Id. M e establish now the corresponding integral formula, 
for which we prove first that if/(r) is any vector function (not 
necessarily continuous) mtegrable along a given curve AB, then 

P 

P being any variable point on that curve, the integral J'/Xr). flf/' 

A 

is a continuous function of the position of P on the curve. 

P 

Denote ^f{r).dfhy P(p). Then if Q is any other point 
A 

Q 

% 

r+€ on the curve, we have j /(;•).,/;— FO +c), and therefore 

A 
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F(y+€) — F(?*) = 



But (see p.l6) 



<Ij 7, where U is the upper limit 


of/.^ for the curve AB and I is the length of the arc PQ. 

Hence | F(rH-€) — F(r) | <U/, 

and therefore | F(r+€)-F(r) | can be made less than any 

V 

arbitrary positive number 8, if only 2 is so chosen that 8 >U/. 



which is always possible because U is supposed to be 


finite. Again, since the curve is supposed to possess a continuous 
tangent at P (p.lO), there is a finite portion of the curve about P 
for which the arc measured from P and the corresponding chord 
increase together. It is possible therefore to take a point P' 

on the euive in such away that the arc PP' <^and also such 

that the arcs corresponding to chords PCI which are less than 
I P P' 1 , are less than the are PP' and less therefore than — . 
It follows that for all vectors c satisfying j c | < | pp' | ^ 

where the are PF< we have i F(r + e)-F(r) | <8, which 

proves our theorem . 

14. We conclude that if/(r) is integrable along any curve 
in a continuous region, and if A is a fixed point and P any 
variable point in the region, then integrating along the various 
curves through A and P, we have any number of functions 

each of which is continuous for a displacement of P on 
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the curve along which the integral is calculated in any ease. 
Under certain conditions however (See §42), of which the 
continuity of/* is one, the intergral is known to be independent 

p 

of the path of integration, and in this case therefore / 

• A 

will define a uni(|ue continuous function F(r) of the position of 
P in space. Assuming these conditions to hold, and assuming 
in particular that /*(?*) is a continuous vector function, we shall 
prove here that /= VF. 

For, if f is continuous, and since in accordance with the 
understanding in § 9, the unit vector t to any curve through 
A and P is supposed to be continuous,/**^ is also a continuous 
function of the position of a j^oint on this curve. Hence, if P JS 
the point r and we consider another point / on the curve in the 
neighbourhood of P, f and t' being the values respectively of 
/‘and t at r', then for any arbitrary 8, a positive number 7 / can 
be found such that 

1/'/-/'^ I <8, if I r'-r I <77, 

which shows that | lies between |/'*^ | +8 and |/*^ I 

— 8 ; every value, in other \vords, oif't in the portion of the 
curve between r and r' lies between \f't\ +8 and \f'^\ 

r+€ 

If € denotes the vector f'—Vy the integral / /’*^^’ hos between 

r 

l[\ ^'•t\ +8] and \f'i | —8], I being the length of the 
arc between r and r . 

>•+€ 

But F(r + e) — F(;*) = / /*■ (h\ 

r 

F(r + €) — F(r) lies betwern \ +8] and I 

But 8 and therefore also I can be taken arbitrarily small; 
and with the arbitrary shortening of I the difference between the 
arc I and the chord | « | becomes arbitrarily small, and the 
direction of c approximates to that of t. We can write therefore 
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F(#* -f c) — F(r) = c \f(r) + 
simultaneous limit zero. 



where S' and 


have 



This result, which is true for all curves tljrousrh A and P and 
true therefore for vectors c drawn in all directions round P 
shows that /= v F. 

If therefore we have any curve in the region, and r,, are 
any two points on the curve, we have 



VF-^/r= J /•^r=F(rJ-F(rO 


IV 

The Linear Vector Function. 


15. The most general vector expression linear in r can 
contain terms only of three possible types, r, aj'b and c-c?', 
n, d, c being’ constant unit vectors. Since r, a.rh and are in 
general non eoplanar, it follows from the theorem of the parallele- 
piped of vectors that the most general linear vector expression 
can be written in the form 


X?'-h /A a.r6 -|- vc X r 

wliere ;a, v are scalar constants. The constants /a, v may 
moreover be incorporated into the constant vectors a and c and 
we write our general linear vector function in the form 

<l> (r)=Xr+a,rh-\-cX7% 
where d only is a unit vector. 

Obviously, ^ (r) is distributive ; 

that is, if} (a-p/?)= 0 (^), 

and further <l> (At) =k (r), whei’e k is any constant. 

16. Theorem , — The surface integi’al of (r) over any closed 
surface S bears a constant ratio to the volume T enclosed by 
the surface, the constant depending only on the function but 
being independent of the particular surface over which we 
integrate. 

To prove this we integrate separately the three terms of 4* W 
over the surface. 


s s 

We know /Ar.d(r=:AJr.do-=3AT. 

To calculate fa.7‘h. dcr we break up the region S into tnin 
cylinders with axes parallel to d. Since the surface is closed, 
each of these cylinders like PQ will have an even number of 
intersections with the surface, as in the usual proof of Green s 
Theorem. It is enough to consider here the ease where there 
are two intersections only, the extension to the general case being 
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obvious as in that proof. If then d(T and da- are the elements of 
surface on S enclosed by the cylinder PQ, we have 

h. da = — h. do-=area of the cross section of the cylinder PQ. 

Let P be the point r, then if -c is the length of the 
cjlinder, Q is the point r+xd, b being a unit vector ; and 
the sum of the contributions of da and d<T to the surface integral 

O 

approximates, as the cross section of the cylinder diminishes, to 

a.(r+a:6) 6. do-' + a. r h, da 
t.e., to a.h jih.da' 

i.e., to a.hdTy where dT is the volume of the cylinder PQ. 
Hence the whole surface integral 

S 

J‘a.r6.dcr=a*6T. 

We may just by the way note from the symmetry of the' 
result that fa,rb,da—^b.ra.da=a’bT!, and this result holds for 
any two arbitrary constant vectors a, .6. That is, for any two 
constant arbitrary vectors a, 6 we have 

a.frh.da=b,la.r d<r=a'6T, which shows, moreover that 
J/*6.do- = 6T and Ja.?*d<r=aT. 

Generally therefore /ra,do’=Ja.?-d<T=aT, a being any 
constant vector. 

To return to our proof now, we have to integrate fc x r.da, 
Putc=ax)3, so that cxr=za.rp—p.ra (p. 6) 

Hence /c x r.da—ia.rP'da^fp,raJa=za'PT-~a’PT=0. 

We have therefoi’e finally 

f0(r).d(r=3XT + a.6T 
^ dor=3X+a.6=D, say, 

which proves our theorem. 

17. The skew surface integral of ^{r) over any closed 

surface S divided by the volume T enclosed by the surface 

is a constant vector, this constant vector depending on the 

function but being independent of the Burface over which we 
perform the integration. 

4 
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Proof. We proved in the last ai’ticle that fcxr.d(r=0, 
being any constant vector. It follo’n-s that fc.r x d<r-0 (p. 6), 

t.e., c.jrxd<T=0, or J’rxda’=0, because c is arbitrary. 

Also, Ja.rh X d(7=b X Sa.rd<x=b X aT [§16] 

Again /(c xr)x d(r~{rc.d<T—Jcr.da- [p. 6J. 

=cT-3cT=-2cT [§16]. 

S S 

Hence /</.(»•) x d<T=:f[\r+a.rb + cx »•] X d<T 

= -(ax6+2c)T=-CT, say ; 


C — — ^ /<^(r) Xd<r= ~ fda-X ^(r) 


being a constant vector, 


our theorem is established. 

D and C which we find here associated with every linear 
vector function, we shall always refer to as the scalar 
vector constants respectively of the linear vector function. 

18. We consider the function now 

cx r 

which is immediately seen to have the same scalar constant 3X+ 
a.b as the original function <(>[)') =Xr-\- a. rh + cx r ; and its vector 
constant is--(a x 6-|-2c) which differs only in sign fi'om the 
vector constant of 

Further, if a, (3 are any two arbitrary vectors, we have 
=a. 4- a'pb + c X y3] 

=/S’[Aa4-6*aa — c X a] if we call the new function 

The two functions <t>(r)=\r + a.rb-\-cXr 

and €l>Xr)—Xr‘\-b'ra~cXr 
may on this account be called conjugate functions. 

With every linear vector function then is associated 

another function characterised by the propery = 

for any two arbitrary vectors and having further 
the same scalar constant as and a vector constant differing 

only in sign' from that of ^(r). 
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19. Since the scalar or vector constant of the sum (or 

the difference) of two linear vector functions is obviously the 

sum (or difference) of the scalar or vector constants of the 

two functions, it follows that the scalar constant of ^(r) + 

is 2D and its vector constant is zero ; and that the 

scalar constant of is zero and its vector constant 

is 2C. 

Obviously again the conjugate of <^(r) +<^>'(0 is itself ; this 

function, that is to say, is self conjugate. And the conjugate 

of <#>(r)-</,'(r)=<^X0-^W = -[<^(r)-<^'(r)], which is the origi- 
nal function with the minus sign prefixed. Such a function 
has been called skew or anti-self-conjugate. 

The function <^(j-)-<^'(r) in full is 

a.7'h — b.7-a-\-2c X r 

= (ax6)xr-f-2cx^- [p. 6) 

= (ax6-f 2c) xr=C xr. 

Hence </>(?•) can be written 

=i[<^>(r) + <#)'(r)] + i[<^(r) - ^^r)] 

=<#>.(r) + i Cxr, where 2<t>,(r) has been written 
for the self conjugate function ^(r)d-<^'(r). 

Any linear vector function ^(r) therefore can be expressed 

the sum of two other functions one of which is self conjugate, 

has the same scalar constant as <<,(r) and no vector constant, 

and the other is skew, has no scalar constant and the same 
vector constant as 

The result ,^(r) =C xr shows moreover that the vector 
constant of all self conjugate functions is zero, 

20. The vector constant of ^.{r) may be exhibited in another 

manner, for which we calculate first the gradient of the scalar 
function . 

Since 8[r.<#,(r)] = (r + 8r).</,(r + Sr) -r.<^(r) 

= (r + 8r) .[<^(r) + ,/,(8r)] — r.<^(r) 

=r.0(8r) + 8r.^(r) + 8r.^(Sr) 

=8r.[<#.'(r) + <#,(r)] + , | Sr | , where r, is 
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a scalar number which has limit zero as | | tends to vanish, 

it follows that V [?’.</>(?*)]=</>(;■) + ^'( 7 '), and we can write 

<t>(r) =1 V[r.^(r)] + X r. 


Integrating now round any plane closed curve, we have 

+ xr.fZr. 

Since ?•.</>(;•) is single valued, the first integral on the right hand 
side is zero, because it is equal to the difference in the values 
of I ?’-</>(^‘) at the same point before and after circuiting, [p.23.] 
Also Jrx^Z?* is twice the vector area enclosed by the curve, 
a fact which becomes obvious by taking a new origin O' in the 
plane of the curve. For if 00'=a, OP=r and OT=p, we have 
r=p+a and dr~dp^ and frxdr=f(p-\-a)xdp. Hence since 
fdp and therefore also fa x dp vanishes, the curve being closed, 
we have p* x dr=fp x dp which is a vector normal to the plane 
of the curve and equal in magnitude to twice its area. 

Thus J* C x rJr^^C.fr x dr=zC.nS, where S 

stands for the area enclosed by the curve and n a unit vector 
along the normal to the plane of the curve. 


The ratio i J’<^(r).^/r=C.w does not then depend on the pai*ti' 

o 


cular curve round which we integrate, but it depends on the 
orientation of the plane of the curve, on the vector n. This 
ratio obviously again attains its maximum value when n is taken 
in the direction of C. The vector constant of <l>{r) then is a 
vector in the direction of the normal to that plane, round any 
curve on which if we calculate the line integral of the ratio 
of this integral to the area of the curve is maximum, and the 
magnitude of the vector constant is equal to this maximn® 
ratio. 

21. There is just one bit of work more in connection with 
linear vector functions before we are ready for the Differential 
Calculus of vector functions. 

If a is any constant vector, a x </»(r) is of course also a linear 
vector function of r. We proceed to find Dj and Cj the scalar 
and vector constants of ax<l>(r). 
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Integrating over any closed surface (enclosing volume T) we 
have by definition 

D,T=JaX </>(r).d<r 
. =a.J<^(r) X d<r 

=— a.C T, Avbere C is the vector constant of 
Dj = — a.C. 

Again, — CiT=J[a x X da 

=J<j>{r)a.da — ja<j>{r).da (p. 6), 

But J‘a<^(r).d<r=aJ<^(r).do-=aDT, D being the scalar constant 
of 

9 

Also f<fi(7')a.da is calculated immediately by breaking up the 
volume into thin cylinders with axes parallel to a, as in §16, 
p. 25, Thus if OP=:y and PQ=.Ta, the sum of the contributions 
of the elements of area da and da' at P and Q approximates, as 
the cross section of the cylinder PQ diminishes, to 

^(r+oja) a.do‘' + <^(r)a.do- 

which again, since ara.do-'s: — xa.da=\o\. of the cylinder 
and *^(r+.ra)=<^(r) + .'C<^(a), 

approximates to <^(a)dT, T being the volume of the cylinder 
PQ. Hence /<^(r)a.d<r=<^(a)T 

Hence finally— CiT=<^(a)T—aDT 

Cj=13a— ^(a). 



The Differential Calculus of Vector Functions. 


2^. The Differential Calculus of the scalar function of a 
•single (scalar) variable concerns itself with the rate of change of 
the function with respect to the variable. In considering in the 
same way the rate of change of a vector function / (?•) of the 
position of a point in space, the first difficulty we meet with 
is that this rate of change is different for the different directions 
in which the point P may be shifted. In fact, the position of 
P being specified in the usual way by the vector r drawn from a 
fixed origin, a change in the position of P of magnitude li 
and in the direction of the unit vector a would be denoted by 
<th, and the change in the value of the function would be 

/(i 4 «/) /()•). The rate of change then in the value of y at 
P for displacement of P in direction a is 

L /(r+a^)— /(}•) . 

h=.o h 

In Gibbs’ notation this is denoted by tf.V/; we shall often 
denote it, — perhaps a little more expressively — also by da/. 

In order that the limit may exist it'is necessary that y should 
be continuous at P in the direction a. For if f is discontinuous 
in this direction, then however small h might be, | y(^‘ -I- 
/(/*) I would be greater than a certain positive number Sand 


therefore 


/(r+aA)— /(?•) 
' h 


can be made greater than any arbi- 


trary positive number, and therefore the limit cannot exist. But 
the continuity of f alone in direction a cannot ensure theexis- 
tence of the limit, for which it is necessary that the fluctuation' 



1 

h 


f{r^ak)^f{r) 


as a function of h should be arbitrarily 


small within a sufficiently small interval on the line a in the 
neighbourhood of P. The continuity therefore is a necessaiy 
though not the sufficient condition for the existence of. the limit 
in question. 
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But in any case where the limit does exist as a definite 
function of a and r, it is clear as in § 10, that that function 

will be linear in a. We may denote therefore L ^y\r-\~ah) — 

fl 

/('■)] by (a, r), or more simply by (a), where 4>{a) stands 
foi a linear vector function of a. The explicit presence of / in 
r) would serve to bring out the fact that the rates of 

change of/(r) are given by different linear vector functions at 
different points P in the region. 

If tb'e limit exists for all directions issuing from the point 
P,— for which it is of course necessary that /‘should be continuous 
at P — we write for any a 

da for fl. V/=<#) (fl), 

where rj is a vector such that 
I I has limit zero as /i tends to vanish. 

Or, since k {a)=<l. (ah) [ § 15 , p. 24 ], if &f denotes the 

vec r increment of / corresponding to the increment Sr of r, 

+ I 8?’ I . 

. 23. If the shift Sr be supposed to take place along a definite 
contmuous curve r=^(0, then we know (§ 4), that /would 
be a continuous function of t along that curve and we would 
wiife from our relation of the last article that 

which, since <#.(8r) is a linear veetor function of 8r, 

p. 24] 

=0 

Further, if specify any two points K and L on the curve 

»■ and/ , and /2 are the values of / at K and L respec- 
• ively, then we have from the Mean value theorem of §6, 


A /i 
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where denotes the value of ^ at some point M ou the 

curve lying between K and L. In other words 


if, of course, a dehnite ^ exists at every point on the curve be- 
tween K and L. 

In particular, if the curve is a straight line in the direction 

of the (unit) vector a and h is the length of KL so that KL=fl^j 
we may write ^(rt, r-\-aOh) 

—hdaf{r+aBh), 

assuming, of course, that dj exists at all points on the line KL. 

24. We shall practically always confine ourselves to func- 
tions which are not only continuous within a certain region, but 
are also such tliat daf or exists at every point P of 

the region for all directions a round that point. If 
construct a sphere of unit radius with P as centre, then 
every point on this sphere will represent a definite direction 
issuing from P, and being known for P would mean that 
corresponding to every point on the unit sphere we know the 
rate of change of / (at P) both in direction and magnitude. 
But as to the rate of change of / at P as a whole, we cannot as 
yet form any definite conception, not at least directly from our 
knowledge of the function <#> at P which only brings before our 
minds a bewildering diversity of the rates of change for the 
infinitely many directions round P. What we naturally do there- 
fore is to have an idea of some sort of average value of ^ (^) 
for these directions average of ^ (a) over the unit sphere 

round P. 

We consider then two kinds of such an average value. 

Since a is a unit vector, the magnitude of the component 
of <l>{a) in direction a is a.^(a). We consider first the average - 
of this magnitude over the unit sphere, which is 

Ja*<j^(a) dS 
S 
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where rfS is the scalar element of area on the surface of the 
sphere at the terminus of the vector a, and S is the whole 
surface S. 


Since the radius of the sphere is unity, S is equal to 4^- and 
is therefore equal to 3T, where T is the volume of the sphere. 
Also adS=(i<T the vector element of area, the normal to the 
sphere at the terminus of a being along a. Hence the average 
value 


where D is the scalar constant of the linear vector function 
^{a). [D is a constant here in the sense of being independent 
of a, but is of course a function or r and is different from point 
to point] . 

The average rate of change of / (r) then in the direction of 
the displacement of P is proportional to the scalar constant of 
and this average rate therefore for each point of the sphere 
may be constructed geometrically by making the sphere bulge 

4 

out uniformly outwards from its centre P by an extra length 

proportional to D. 


To have an idea now of the average value of the tangential 
component of </>(a) for all the points of the sphere, we consider 
naturally the average value over the sphere of the moment of 
about the centre P. This moment being a x our 
average value 


JaX<li(a)d8 /do-x^(a) 

S S 



where C is the vector constant of in the present instance 

of course the constant being a function of r* The average 
moment therefore is in the direction of C and in magnitude is 
proportional to that of C. It follows that the average tangential 
component of 4>{a) on the sphere, — the average of the component 
that is to say, perpendicular to a,— ^is perpendicular to C and 
in magnitude is proportional to that of C. The vector constant 
of ^(a) affords us a knowledge of the average rate of change 
5 
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y*(^) for any small displacement of r, perpendicular to that 
displacement. 


The scalar and vector constants of <l>(a) therefore may be 
regarded as supplying us with a basis of comparison of the rates 
of change of /’(r) at the various points of the region, and serve 

in a sense the same purpose that is served by our old ^ in the 


case of scalar function of a single (scalar) variable ,)*. 

25. The scalar D and the vector C being then so fundamen- 
tal in the Differential Calculus of Vector Functions, we hasten 
to exhibit them directly in terms of the function y(r) to which 
they belong, and we shall find incidentally how they are ulti- 
mately identified with the well known Diversfenee and Curl of 

o 

the vector function /(r), 

26. Divergence, We integrate //. over any small closed 
surface surrounding the point P. At any point Q on this 

surface, PQ being Sr, the value of/ 


where ^ stands for | 8r | and ?? is a vector such that | v * 
zero limit as h approaches zero. 

Hence //.6?o'=J[/p+^ (Sr) 

But Ifp.dtT =/p.j‘(^cr=0, the surface being closed ; 

f<j>(Sr).dtr='DT, D being the scalar constant of <t> (M 
regarded as a linear vector function of Sr and T the volume 
enclosed by the closed surface. 

I rj I being the scalar magnitude of dtr 


< rj fhds, where yj is the greatest value of | V 


$ 


But fAdS=KT, where k is a.finite number 
Therefore, /yiiy.do* <^'cT. 


Hence — //.dir— D < t/k. 
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Now let the surface shrink up to the point P in any manner. 
Then since when h approaches zero, all | >? | ’s and therefore 

also Tf tend to limit zero, it follows that L exists and is 

equal to D. We have the following definition then — 

Enclose the point P by any small closed surface and calculate 


the integral J /.do-. If the limit of — / /.do* exists as the surface 

shrinks up to the point P, which limit moreover is independent 
of the original surface and of the manner of its approach to 
zero, then this limit is called the divergence of /(r) at P. This 
limit exists if <^(a) exists at P, and in this ease the divergence 
of /(r) is equal to the scalar constant of </>(«), and may there- 
fore be taken (but for the constant factor i) as the measure of 
the average rate of change of f(r) corresponding to any dis- 
placement of P in the direction of that displacement, the 
average being taken for all directions round P. 


We shall always denote the divergence of y{r) by 

27, It is perhaps possible for the divergence of a given 
function to exist at a given point P without necessarily 

existing there. Directly, the necessary and sufficient condition 
for the existence of the divergence at P is that within a suffi- 

g 

ciently small neighbourhood of P, the function ^S/.dtr should 


vary continuously for continuous variations of the surface S. 
In other words, given any arbitrary positive num ber 8, it should 
be possible to find a positive number such that S,, S, being 
any two closed surfaces round P and T^, T, the volumes 

enclosed by them, 


1 S. 1 s, 

1 s 

should be less than 8, whenever the surfaces S^, S, are entirely 
contained within the sphere with centre P and radius rj. It is 
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a matter for investigation now how far this condition necessarily 
implies the existence of </>(«) at P. We consider however just 
now only those functions for which exists at every point 
and for which therefore there is no question as to the existence 
of the divergence. 

28. Cxirl, — Integrating J'/x(/<r over any closed surface in 

the same way^ we have 

//x do-=:{[fp-^<j>(Br) + hr/] X d<r. 

//p X t?o-=/p x/iio-=0, the surface being closed ; 

j <t>(8r) X d(r= — CT, C being the vector constant of <t>(8^') 
as a linear vector function of Sr. 

Also fhy)Xd<x'j^fh\T)\dS 

that is,< 


^ //xdcr-f-C< ??K, and therefore in the limit when the 


surface shrinks up to the point P, we have 

Sfxda=-C. 


\ In general, if we find that the limit oi d<rxf exists as 

the surface shrinks up to the point P, the limit so obtained 
being independent of the original surface and of the sequence 
of forms taken by it during its approach to zero, then this limit 
is called the curl of the function /(r) at P. What we have 
proved above shows therefore that if <f>(a) exists at P, the curl 
does so too and in this ease is equal to the vector constant of 
4>(a), For any displacement of P then, the average (for dis- 
placements in all directions round P) of the component rate of 
change of /(r) perpendicular to the direction of the displace- 
ment is perpendicular to the curl, and the magnitude of the 
average is proportional to that of the curl. 
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29. The curl ma}' also be exhibited in another (the more 
usual) manner corresponding to the property of the vector 
constant indicated in § 20, p. 28. 

Thus, integrating J/. dr round any small closed plane curve 
surrounding P, we have 

J/.dt- = /[/p+<;>(8r) + ^].d7-. 

But J/p'd)'=/p*J'dr=0, the curve being closed. 

‘dr=C -nS, S being the area enclosed by the curve 
and n a unit vector normal to the curve. [§ 20]. 

Also, Ihri'd/r > J/i 1 7} I ds^ds being the magnitude of the vector 
element of arc dr ; 

and fh \ yj \ ds <'^ J?idS, being as before the greatest | t? | , 
and also J/^d^=^cS, where k is a finite number. 

lhyfdr< ^kS, 


and we have i J/'dr— C*n < « k. 

S ' 


Hence, as the curve contracts to a point in any manner, 


1 


g J /*rfr approaches the unique limit C*/z. The limit moreover 


exists for all aspects of the plane area, for all vectors ?t. These 
various limits [for the different aspects of the plane have again a 
maximum value when n is in the .direction of €, that is, when 
the plane is taken perpendicular to C and the magnitude of this 
maximum value is equal to that of C. 

In case therefore <l>(a) exists for the function /* (?) at a point 
^ » it is indifferent whether we define the curl as we have already 
done it, or use the following definition. — 

Having described any plane closed curve surrounding the 

point P, if we find that the limit of ^ as the curve con- 

s 

tracts to the point P, exists and is independent of the form of the 
curve and of the manner of its approach to zero, but dependent 
on the orientation of the plane of the curve, and if further, as 
this orientation is varied, the various similar limits so obtained 
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for the different orientations all exist and acquire a maximum 
value for a certain orientation, then a vector drawn perpendi- 
cular to the particular plane which gives us the maximum value 
and equal in magnitude to this maximum value is the curl of 
/(r) at P. 

We shall denote the curl of f{r) by v x f{r) 

30. Some Transformation Formulae , — The explicit recogni- 
tion of the divergence and curl of a given function as the scalar 
and vector constants respectively of the linear vector function 
^ Oj which defines the rate of change of the function for 
displacement in any direction a, considerably facilitates the mani- 
pulation of these operators in practical work. This is what we 
proceed to illustrate. 

We shall in this article denote by n a continuous scalar 
function possessing a gradient at every point within the region 
considered, and U and F will stand for two continuous vector 
functions of which the rates of change at any point P will be 
denoted by the linear vector functions ip (a, r) and ^ {a,r) res- 
pectively, so that the scalar and vector constants of will 
be and x U respectively, and those of tt>{a) will beV 
and V X V respectively. 


(^) To show now that 

V.(wV)=« V.V -h Vw.V. 
and V X (mV)=7*V x V-h Vwx V, 


Proof. Q being a point in the neighbourhood of P(PQ=®^)> 
if the values of u and V at Q are u+Bu and V-f-SV, the rate of 
change of wV for a displacement in direction a 


= L i [(w + Sw) (V-f-SV) — ttV] 
h=o ^ 



[w8V-f- VS m + SttSV] 


= -f Va. V w 

Hut V * (wV) and V x (uY) are the scalar and vector constants of 
this rate regaaxled as a linear vector function of a. Remembering 
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therefore that the scalar constant of a,rh is a.b [p. 25] and that its 
vector constant is ax 6 [p. 26], we have immediately 

V(wV)=tiVV+Vtt-V 

and Vx(wV)=mV xV+VwxV. 

(m) To show that 

V(UxV)=V. VxU-U- VxV 

and V X (Ux V)=UV ■ V-W •U+^(V)-</»(U) 

Proof, The rate of change of U x V for the displacement ah 
of the point P 

% 

=Li[(U+8U)x(V+SV)-UxV] 

■ =Li [Ux8V+8UxV+8Ux8V] 

=Ux<^(a)+i/^(o)xV; 

and we have to find the scalar and vector constants of this as a 
linear vector function of o. 

We lecall [§ 21, p, 29] that the scalar and vector .constants of 
are a-C and Da — <^(a) respectively. Hence 

V tUxV) = -U- VxV+V. VxU 
*nd Vx(UxV)=Uv-V-d)(U)-VvU+i/r(V). 

=uvv-u- vv-w-u+v- VU. 

(»»t) We may quite easily prove also the formula for V(U-)V 
given in Gibbs, Vector Analysis, p. 157, viz. 

V(UV}=Ux(VxV)+Vx(VxU)+U- VV + V- VU. 

Thus the rate of change of U • V for the displacement aA of the 
point P 

[(U+SU)-(V+«V)-U V] 
=I^^tU-8V+8U.V+SU-8V] 

=U-.^(a)+V-^(o). 
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Now we know that if this rate of change can be written in the 
form a.G., then G=V(U V.) 

But U' (^(a)+V i/^(a) 

=U- [<^(a)— ^'(a)] + v- [./.(a)-i/^'(a)] + uy(a) + vy(0) 

=U- (VxV)xa+V- (VxU)xa+o-<^(U)+a-i/^(V) [p. 27] 
=a- [Ux (V X V) + Vx(V xU) + <^(U)+^(V)] 

Hence V(U ■ V)=U x ( V x V)+V x ( V x U) + ^(U)+i/'(V) 

=Ux(V xV) + Vx(VxU) + U'VV+V-VU 

We could write the same result in a more compact form. 

Since U'<^(o) + V'i/'(a) could be written directly 

X a-ifr'{Y), we have 
V(U-V)=<^'(U) + f(V). 

In particular, V((i V)=<^'(“). « ^ constant vector. 


A short note on Bilinear Vector functions. 

31. Before passing on to the Second Derivatives of the 
vector function /(r), it would be necessary to consider very 
briefly what are called the Bilinear Vector Functions. 

A vector function of two variable vectors linear in both 


called a bilinear vector function. 

Generaljy, a vector function of n variable vectors linear 

in all of them is called an w-linear vector function. 

A bilinear vector function is said to be symmetrical if it 
remains the same when the two vectors are interchange . 
r r denote the two variable vectors, the general symme ri 

bilinear vector function can contain only terms of the W 

where X is a constant vector and is a self^o ] g 

linear vector function, so that e 

therefore for the symmetrical bilinear vector function 


+r'/',(r')X, 

or, = 

where all the functions i/' are self conjugate. 
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32, Since the scalar constant of a'rb is a b [p. 25], the 
scalar constant of regarded as a linear vector function 

of r alone, is or since the functions i/* are all self- 

conjugate, this scalar constant. 


of which again, regarded as a function of ?•', the gi-adient is 

S-ACA). 

It 18 now obvious a priori from the symmetry of {, ■,)■') and 

it is verified immediately also-ihat if we had calculated the 
scalar constant of regarded as a function of >■' and then 

had obtained the gradient of this scalar constant with respect 
to n we would have got the same result Hence, without 

ambiguity, we may refer to as the gradient of the scalar 

constant of </>(r,r'). We shall denote §i/r(X) by f 

33. Since again the vector constant of a.rb is ax 6 [p. 26) 

Uie ^^ector constant of <^(r,,-'), regarded as a function of r, is 

X A. We want to write down now the scalar and vector 
constants of this vector constant 5^(r')xA regarded as 

a function of r . We recall that the scalai- and vector constants 

o ax<t>^r) are a'C and Da— <^(a) respectively, where D and C 

are the scalar and vector constants respectively of d>(r) Hence 
since the functions ^(r') are all self conjugate and therefore’ 

ot X A as a function of r' is zero. 


And its vector constant 


-['Ai(A,)+^,(A,)+ ]-[A,Dj+A,D,+ 


]. 


' 

^ D,...are the scalar constants respectively of 

That 18, the required vector constant 
=5WA)— SAD=f- — gXD 

tha aame if we 

oW "®gardeda« a func 

tht vlt ’ * constants of 

this vector constant as a function of r, 

6 
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We inay say then without ambiguity that the scalar constant 
of the vector constant of Is zero and that the vector con- 
stant of the vector constant of 0(r,/) = f - We shall denote 

this last by C'. 

34. For the bilinear vector function we have -now to consider 
two conjugate functions. If we regard it as a function of r alone 
we have one conjugate function, and we have another when we 
regard it as a function of alone. We denote these two conju- 
gates by <!> and <f>'r’(ry) respectively. These conjugate 

functions are not necessarily symmetrical. Remembering that 
the conjugate of a'rb is 6*ra, we have in fact 

if 5r'*i/r(r)X 

and — 

and we propose now to seek for constants like f and C' from 
these conjugate functions. 

Since a linear vector function (of one vector) and its 
conjugate have the same scalar constant and vector constants 
difFering only in sign, it is obvious that the scalar and vector 
constants of <j>\ regarded as a function of r and of <!>'/ regarded 
as a function of r' could be written down immediately from 
the results we have already worked out for (r, ?*'), but they 
would not obviously also furnish us with anything new. Also 

is only with ;• and r interchanged. We have to cal- 
culate therefore only the scalar and vector constants of <t>'r 
regarded as a function of r'. 

Since the functions are all self conjugate, the vector 
constant in question is zero immediately. And the scalar 
constant = 5X*rD = r*^XD, the gradient of which with respect 
to r is ^XD. We denote ^XD by f Thus f '=5XD is gra- 
dient with respect to r of the scalar constant of regarded 
as a function of r' and is also the gradient with respect to 
r' of the scalar constant of regarded as a function of r. 

We thus have two independent constants of the symmetrical 
bilinear vector function, viz, p and p ' and a third C' dedneible 
from them. 
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If tbe function is 

r=5^(X) 

r'=5XD 


and C'=r-r' 


86. Now we are in a position to consider repeated opera- 
tions of the derivative operators. 

Denoting L the rate of change of f(r) 

k = Q 

at P for a small displacement of P in direction a by (a, r), 

we consider first of all the mte of change of this function '<l>,{a,r) 

for any displacement of P. If a'h' is this new displacement, 

a being a unit vector and // a small positive number, then the 
rate of change of (a, r) 


/ i '=0 


=L 
^' = 0 



= L 


]; L \{AT + a'h'^ah)-f{r-\-a'h')-f{r+ah)^f{r)-]. 


A'=0 > 1=0 

Assuming that a unique Umit exists as h, h' approach zero 

we conclude as in §10, p. 17 that this limit is a vector function 

linear both in a and a'. We denote this bilinear vector function 

by ^8(a', a, r) and call it the second differential linear vector 
function for/(r),-^i («, r) being the Erst. 
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By definition then 


In the same 


way, 


{a, r), 

Wt r). 


d,dJf=L-^ 

h 


1 


[/('■ + + «/^) -/(yH- */(; H- 4^^^ 


/^ = 0 /i' = 0 

and this would be denoted by <l>^ (a, a\ r). 

We see that djdaf and d^djf differ only in the order in 
which h and 1i are made to approach zero and under certain 
circumstances, which may be investigated, the limit operations 
are commutative and then we should have 


djd^f—dadj/. 

We prove here only that in case daf, djdaf d^d^f 
all exist and are continuous within a finite region round P, the 
commutative property of the limit operations in question 
certainly holds and we have dadaf—dad„'f 


For, applying the Mean VaJue Theorem of §23, P. 32 to 
the function 


f (r). 


we have [ f{r-\-a'Ii-\-ah)~-f{r-\-ah)'] — /(r)] 

= hd„\_f{r-\-ah'-\-a6?i)-~f{r-\-a6h)],0 being a positive 
proper fi-action, 

=hdJidJ\_f{r-\-a'd'h' + a0h']y applying the same Mean 
Value Theorem to f{T-\-a6h)i 0' being some other positive proper 
fraction. 


In the same way we have 

[/(r+a'>^' + a>5)-/(r + ^7/)] - [/{r&ak) 
^h'd.'lf{r-\-ah+a'e\h')-f{r+a'e\^ //')] 

m 

being also positive proper fractions. 
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Provided only then r-^ahy and r ah -{■ a* k' lie within 

the region where our conditions hold, we have it that whatever 
hf h! might be, 

f{r -f aJh! + ah) -/(r + a'h’) — /’(r + ah) +f{r) 

= // V/„ [ /(r H- // + 1 /O ] ■ 


Since again d J a' & h' adh) and d J d J[r a! e\h' ^ 
ad^h) are supposed to be continuous, they approach the same 
limit at r, and we have djdaf=:dadaf\ or =^ 2 (^)^^ 0 • 

In other words, is a symmetrical bilinear vector function of 

j 


a, a . 


36, Since by definition regarded as a linear 

vector function of a gives us the rates of change of </>i (a,?*) 
for the directions a\ the divergence and curl of this latter 
function are the scalar and vector constants respectively of 
<#>,(» as a linear function of a'. We proceed to show now 
how the second derived functions of /(r) can be obtained from 
just as we obtained our first derived functions — the diver- 
gence and curl of /(r)— from Since the divergence of f{r) 

is a scalar, we can have its gradient, and the curl being a vector, 

we can have its divergence and curl. We consider these in 

order. 


37. Let D stand for the divergence of f(r) ; D=:V-/- If Q 

is a point in the neighbourhood of P, such that PQ=:aA, we have 
for the divergence of /'at Q, 

D+SD=V.(/+80 = V./+V.8/ 

/. 3D = V.8/= (a/t,r) 

a,VD=Liv 

fe=o ^ 

0 being a 

linear vector function of a, 



46 


VECTOR CALCULUS 


But is the scalar constant of regarded as a 

linear function of a' and is therefore equal to a* (" [§ 32], Hence 
a.VD=a.r , and since a is arbitrary, we have 

vD=r, 

where p is the gi'adient of the scalar constant of 

38. In the same way, if C is the curl of / at P and CH-5C is 
the curl at Q, PQ as before being ah^ we have 

C=Vx/, 

C+8C = Vx(/+8/) 

8C = V x8/=V X [<^,(aA,r)+/iT 7 ] 

and daC = L-^V X [^j(a^,r) + 4?y] 

h=o " 

= V 

which we know is the vector constant of (a',ff,r) regarded as 
a linear vector function of a'. 

Hence V.C and V xC, which are the scalar and vector con- 
stants respectively of d„C as a linear function of a, are respectively 
the scalar and vector constants of the vector constant of 
The former we have px’oved to be zero [§ 33, p. 41], 

that is, V. V x/ is always zero ; 

and the latter was found to be f ~ T ^ i® 

Vx(Vx/) = VD-r', 

where j" ' is the gradient with respect to a of the scalar constant 
of regarded as a funtion ofa’» We shall presently 

find an interpretation of f ' directly in terms of /. 

39. It was seen [§ 30, (iii), p. 40] that if 

then where <l>i'(a,r) is the linear vector 

function of a conjugate to ^i(a,r). 

Hence V.V(a* /)= 

that is, V-V(«-/) is the scalar constant of regard- 

ed as a function of a'. 
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Now <;,'<^/(a,r) = Li[<^/(a,)- + a'V)-,^,'(a,r)] 
Therefore, being an arbitrary constant vector, 

{a, r)=h~[P 4^' (a, r+ a' h')-p.,j, a, r)] 

h'=o ^ 

(Aj- 4- a'V) — a.<^j t ^,r)] 


(A^*+ a'V) — (/3,r)] 

fl 0 

=a.<^,(o',/5,r). 

Hence rf.'<^,'(a,r) is the conjugate of ’i>^{a'a,r) with respect to 
a- That is, 

„ (a',a,r) 

Hence V. V(a./)=a- f '. [§ 34, p. 42] 

That 18 , the divergence of the gradient of the (scalar) magni- 
tude of the component of / in any direction is equal to the com- 
ponent of (" ' in that direction. 


40. Generally, starting from any scalar function 

have two second derivatives, the devergence and curl 

ent G= Vw. 


«, we ma}’ 
of its gradi- 


If G + SG devote the gradient at Q, PQ=aA, we have 

G+8G=V(«+8tt) = Vw+ v(Stt) 

8G = V (8m) = V [^ • Vm + Kr{] 

= V [Aa.G + hri] 

/. d„G=Liv[Aa.G+A?;]=zV(aG). 

The divergence and curl of G now are respectively the scalar 
and vector constant, of d.G as a linear vector function of a 

For shoi^ess- sake the divergence of G, that is, v (Val 

IS always devoted by V«. We prove now that the ctT of 
G IB always zero, 
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For denoting c?„G by <t>(a) for a moment, we know that 
V(flt.G)=</»'(a) ; and we have just shown that V(a.G) 

In other words, <^(a) is a self conjugate function and therefore 
its vector constant is zero. 


Hence, V x ( Vw)=0. 

♦ 

41. We have thus obtained two second derivatives of a scalar 
function ?/, viz., V X ( V«) and V'(V«) or V*w, of which the first 
vanishes for all functions u. 


For the vector function /(r) we got three second derivatives 
V( V'/)» V X (V x/) and V'( V x/), of which the last vanishes 
for all functions / and the other two correspond to two of the 
invariants f and C' of the symmetrical bilinear vector function 
a, r). 


But there was a third invariant p ' of (a\ a, r) which we 
saw was related to f by the relation 


VV(a*/)=ar'or V*(a-/)=aT'» 


showing that V®(®’/) Is maximum when a is taken in the direction 
of r ' magnitude of this maximum value is equal to the 

tensor of f This suggests the following definition of a fourth 
second derivative of the vector function /(r). — It is a vector of 
which the direction is that along which if we calculate the 
component of/(r), the divei'gence of the gradient of the (scalar) 
magnitude of this component is maximum, and of which the 
magnitude is this maximum value. This, as the relation 
y*(ay)=aT ' verifies immediately, is of course that old deri- 
vative which is so familiar to us in its Cartesian form 






or tV*w4-iV*v+AV*«’, 

«, V, w being the Cartesian components of/(r) and as. usual 

unit vectors along the axes. But it is certainly significant how,. 
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without consciously seeking for it, we arrive at it all the same 
from an unbiassed and straightforward examination of 

a, r). 

We obviousl/ require now a new notation for this second 
derivative, for it is not dedueible by any repetition of the first 
derivative operators, vt:., the gradient, divergence and curl. As 
all writers on V’ector Analysis, — not excluding mathematicians 
like Silberstein (Vectorial Mechanics, Chap. I) who constantly 
advocate the exclusion of Cartesian decomposition from Vector 
Analysis, — have always defined it by its Cartesian expression 
which immediately suggest for it the notation V */(?*), this 
V*/(r) is the notation that is invariably employed. With this 
notation now we write 

V’/=r': 

V’‘(a-/)=a-V‘/ 

and further, V x ( V x/) = V ( V/) — V */• 


7 



VI. 


Integration Theorems. 


42. The eharaeteristie properties of the divergence and curl 
lead almost immediately to the integration theorems 

S T 

S T 

(n) ffxdo'=^ /V Xfdr , 

the integi-ations extending over the surface S and through the 
volume T of any finite closed region ; and again 

(m) { fJr={V xf.d<T , 

the surface integral here extending over the surface of any 
finite unclosed region and the line integral round the contour of 
the unclosed surface. The function /{r) in all cases is supposed 
to be finite, single-valued and continuous at all points of the 
region of integration. 

To prove the first theorem : — f'or any sub-division of the 
region into smaller closed volumes, we know by the usual 
argument of the integrals cancelling over the interfaces, (the 
continuity of/ ensuring the equality of the values of / at cor- 
responding points on the two sides of an interface) that 


S 

J f,d<r= 


r=n Br 
r=l 


S,. denoting the whole surface of any one of the sub-regions. 

If now we have a convergent system of sub-divisions such 
that at any stage the greatest diameters of all the sub-regions are 
less than any arbitrary number then 

S. 

f f.d(r = DrTr-^Vrrr f § ^6, ] 

where D r is the divergence of / at some point within S,, ,Tr 
the volume enclosed by S,. and is a number having limit *ero 
as h tends to vanish. Hence 
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S 

Let us pass to the limit now as the sub-division advances and 
h therefore diminishes indefinitely. 

Now is alwajs less than i) rj being the greatest 

fr < >? T. 

But T is supposed to be finite and in the limit r) =0. 
in the limit t, =0. 

T T 

Also, by definition, the limit of t,. is / Ddr or JV./dr. 

S T 

Hence 1/^0-= 

In precisely the same way we prove that 

S T 

//xd(r=- /Vx/(ir. 

lo prove the third theorem, we break Up, as usual, the un- 

closed surface by a network of closed cui-ves. Then having a 
definite convention as to the sense in which the line integrals are 
to be calculated round these curves, we prove first of all in the 
usual way that the line integral round the original contour is 
equal to the sum of the line integi-als round the closed curves 
that have been drawn on the surface. Hence, using the same 
sort of argument as used above for proving the first theorem, 
and by a reference to §29, we prove quite easily that 

lf.dr=l V y.f.d<T. 

Corollary 1. For a closed surface /V xf.d<T=0, and this 
affords another proof of the theorem V.V x/=0 [§ 38, p. 46] ; 
for the closed surface may be taken as small as we please, and 
then integration theorem (t) will prove the result. 

Corollary 2. The line integral round every closed cui-ve in 
the region will vanish, if and only if V x/ is zero at every point 
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of the region. But if the line integi’al round any closed curve 
like AD PEA drawn through the two points A and P is zero, 

P 

that means that the line integrals J fAr are the same, whether 

A 

we use the path ADP or the path AEP. Similarly if we 
connect A and P by any other path like AFP, then since 
the line integi'al round the closed curve ADPFA is zero also, 
it follows that the line integral for the path AFP is equal 
to that for the path ADP. We conclude that the condition 

P 

(promised in §14, p; ^2) that If dr may be independent of the 

A 

path of integration and may therefore define a unique function 
of P, (A being a fixed point) is that the curl of f should be zero 
at every point of the region considered. 

43. If we in the theorem (i) of the last article, 

a being any constant vector and u a continuous scalar function 
possessing a gradient at every point of the region considered, 
we have 

/ ua.d<r=f V-(®w)dr 

=/ a.'^udr ; [§ 30, (i) p. 38] 

that is, for any arbitrary constant vector we have 
a.ftid<r = a./ V wdr. 

It follows that /wd(r=/V«dT. 

Again putting/'=aw- in the theorem (nt) of the last article, 
we have 

/« a. dr=/V X (ua). da 

=:/VwXa, do- [§ 30, (i) ] 

=/a. d<rX Vw, 

the integrals extending round the contour and over the surface 
respectively of any unclosed surface. 

Since a is a constant vector, we may write 

a. /wdr=a. jdaX V« 

/, lud^r—ldaX V«, « being arbitrary. 
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44. Differential of an integral . — Since the calculation of 
the gradient, divergence and curl all depend on that of the 
differential of the function considered, it is necessary to have 
formulae for the differentials of functions given in the form of 
integrals, before we can get the result of operation on them by 
these derivative operators. 

Let / (r, X) be a vector or scalar function, X being an 
arbitrary (vector) parameter of the function. 

Consider the volume integral // (r, X) dr of the function, 
the region of integi-ation being bounded by the surface 
P(r, X)=0, where F is a scalar function. The inteeral of course 
is a function of X alone, say x{X). 

Imagine now a small increment 8x to be given to X and 
let T and t' denote the volumes bounded by P(/*, X) = 0 and 
X+8X) = 0 respectively. Then 

/(’•»X + 8X)dr-/’’ /(r, X)rfT 

=f fir, X + 8X)dT+f fir, X+SX) rfr-f /(r,X) dr 
-f [/(nX+SX)-/(r, X)]i^T+/^ /(r,X+SX)dT, 

/ 

T 

/('■> denoting that this integral is to be taken in 

the region between the surfaces P(r, \) = 0 and F(;-, \+8 a) = 0, 

d T being reckoned positive or negative according as it is outside 
or inside the surface F(?', X)=0. 

Now as I 8X I becomes smaller, the surface F(r, X+8x) = 0 
approaches F(r,X)=:0 and the volume between them tends to 
become a thin shell distributed over this last surface; and 
if r and r+8r are corresponding points on the surfaces 
F(rA)=sO and F(v^ -r SA.) = 0, the volume dr of this shell resting 
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on tbe element of area dtr at r on the surface F(;',X)=u 
tends to 

hr.d<T or ^rjtdS 
that is, to dS. 

I VrF I 


M the unit vector alon^ the normal to F(r,K) = 0 being ^ . , 

^1 VrF I 

if VrF denotes the gradient of F(y, X) regarded as a function 
of r alone. 

Hence, as | 8A | diminishes, 

/(r, \+S\) dr 

g 

tends to / / (r, \+8X) dS, the surface integral being 

I VfJ^ I 

taken over the surface S of F(r, X) =0 

To express this integral now explicitly in terms of Sx, we note 
genei*ally that 

F(r + 5r, X + SX)-F(r, X) 

=F(r+^r, X+8X)-F(r, X+SX) + F(r, X+SX)-F0*,X) 

= Sr.Vr F(r+^5r, X + SX)-fSX. VAF(r, X+m) [§ 12, p. 20.] 

which, by sufficiently diminishing ( Bk | and | Sr ( may be 
made to approximate, to any arbitrary degree of accuracy, to 

Sr. Vr F(r, X) + 8 X.Va F(r, X), 

if of course both VrF and VaF are supposed to be continuous 
functions of r and X, 

If now r+Sr be supposed to be a point on the surf^e 
X+SX, in the neighbourhood of the point r on the surface X, 
then 

F(r+5r, X+SX)=0 and F(r, X)=0, 

and therefore Sr. VrF (r,X)+8X.VxF(r,X) can be made arbitrarily 
small. 
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Hence the integral / / (r, X + 8 X)(?t further approximates to 


V ,r 


the degree of approximation remaining the same. 
We write finally therefore 

f{r, k)dT 

= f [/(^X+8X)-/(r,X)]dT 

g 

- / /(.,X+8\) dS, 


our assumptions about the function F being that both v r and 
VkF are continuous functions in each of the vectors r and'x. 

-tS. Suppose now / to be a continuous scalar function 
possessing a gradient at every point within the region of inte- 
gration. Then by the Mean Value Theorem of § 12, p. 20, 

T T 

«/ f{r,\)dr=S Vx/(r,\-(.e.SX)]rfT 
S 


-/ f/(a)+8x. Ax/(r,x-t.e,8x)]^^^^i as . 

being a positive proper fraction. 

If we further assume the continuity of Vx / in X, then since 
the surface S and the volume T are both supposed to be finite, 
the difference between the right hand side of the last equation and 

T. S 

S Vx/(r, X;]dT-//(r, X) VxF SX ^g 

I VrP I ’ 

that is, the difference between this right hand side and 

Vx/<iT-// dS] 

I VrP j ■' 

will have limit zero as | 8X | tends to vanish. Hence 

T S 

Vx/ f((r, X)dT=/ Vx/*--// -7^ Jg 

' I VrF I 
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46. Suppose next that f is a continuous vector function. 
Let daf denote the i*ate of change of /*, regarded as a function 
of X, for an increment of A in the direction of the unit 
vector a. 


If we assume now that d^f \% a continuous function of X, 
then using the Mean Value Theorem of § 23, p. 32, and 
remembering that the surface S and voIuraeT are finite, we 
have, by an arscnment similar to that of the last article. 


daXrX)=/c/a/rfr-//^^^^^iS. 

But the divergence and curl of xW are the scalar and 
vector constants respectively of daxQ^^ regarded as a linear 
vector function of a. 

Now since the scalar constant of the sum of any number of 
linear functions is the sum of the scalar constants of those 
functions, and since further the scalar constant of da.f is 
V ■/, we have the scalar constant of jdafdT =/Va.'/J^T • Hence, 
[remembering that the scalar constant of a-rb=a'b (p. 25)]. 

V •x(A) = Va •//('•, X)rfT=/Vx./rfT-/ 

Similarly, 

AxX/ /(r,X)dr=/ VaX/^^t-/ 

47. If X is not involved in the equation of the bounding 
surface, and its variation therefore does not affect the region of 
integration, we have simply, for a scalar function ^/(r,X) 

V;, 

and for a vector function y^(r,X) 

V\‘ 

and VaX^ fd/r = j 

the restrictions imposed on / being the same, viz., that both/ 
and ^should be continuous functions. 
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48. We conclude our vector calculus here. A much greater 
elaboration of details would certainly have been necessary, if it 
had been intended to present the subject with any degree of 
completeness for purposes of practical applications. There are 
certain obvious extensions moreover which suggest themselves 
immediately from the work done here in the foregoino- pawes • 

^ o r o ^ 

tor example, 

(0 The discussion of the trilinear, and generally of the 

«-linear symmetrical vector fuetions with a view to discovering 

the higher derivative operators. 

(m) The discussion of the improper integrals,— in particular 
integrals of functions having one or more infinities in (he region 
of integration, and of finite functions through regions extending 

to infinity, and the differentials of such integrals ; and of Poisson’s 
equation for vector functions. 

(m) The consideration of what Gibbs calls the determinant 
of the linear vector function, and the development therefrom of 
the Jacobian and Hessian of vector functions on the same lines 
that have been adopted here for developing the ideas of divergence 

and curl from what we have called the scalar and vector eon- 
slants of the linear vector function. 

But the object throughout the present paper has been to 
bring into as great a prominence as possible the one idea that 
the concepts of divergence and curl of a vector function,— which 
we are always in the habit of thinking of in terms of those 
physical ideas that gave rise to them,— do also form, quite apart 
from their physical interpretations, the fundamental notions of 
the Abstract Calculus of Vectors, and supply us with a counter- 
part of the differential co-efl5cient of a scalar function in the very 
real sense of giving us a basis of comparison of the rates of 
change of the vector function from point to point of the field ; 
and further that this new mode of viewing them introduces eon-’ 

siderable simplicity in the practical work of manipulation of 
these operators. 



PART II 


The Steady Motion oe a Solid under no Forces 

IN Liquid extending to Infinity. 

An attempt is made in this Part II to apply vector 
methods to the above problem. It is just likely that it will 
be ^ found to contain, especially towards the end, some new 
results which have not yet been worked out either with the 
cartesian or with vector calculus. 

1. The origin O being fixed in the solid, if we denote by 
the vectors R and G the force and couple constituents of the 
“impulse” that would, at any instant, produce from rest the 
motion of the system consisting of the solid and the liquid, 
and by the vectors V and W the linear and angular velocity 
components of the solid, then arguing as in LamFs Hydro- 
dynamics, Chap. VI, §§ 119, 120, we have for the equations 
of the motion of system. 


I 


dR 

(It 


-RxW= 




RxV-Gx W 



0 ) 


where X represent the force and couple constituents of the 
extraneous forces, the left hand sides being the rates of change 
of R, G when the 'origin system’ [see Silberstein’s Vectorial 
Mechanics, foot note, p. 69] has the velocities V and W. 

If T is the kinetic energy of the system, we have 


and further 


2T = V.R + W.G I 

R = V.Tand G=V«T, ) 


... («) 


where Vt,T and V»T denote the gradients cf T regarded as a 
function of V alone and W alone respectively. [Lamb § 122.] 
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2. We have to express first T in terms of V and W and 
then equations (2) will give us expressions for R and G in terms 
of V and W. 

If we put Tj for the kinetic energy of the liquid motion 
alone and Tg for the kinetic energy of the solid, T = Tj +T 2 and 
we calculate Tj and Tg separately. 

3. To calculate Tj — If U is the velocity potential of the 

liquid motion, we have, if we take the density of the liquid 
to be unity. 

2T,=JUvU.^/a 

where the integration extends over the surface of the moving 
solid. 

Now the velocity potential U satisfies the following 
conditions 

(i) V^U=0 at all points of the liquid ; 

(ii) M.vU=«. (V+Wx?) at any point P on the surface 

of the solid, r denoting the vector OP, and n being a unit vector 
along the outward normal to the surface at P. —For the velocity 
of the liquid at the same point is VU, and the normal 
components of these two velocities must be the same. 

(iii) vU=0 at infinity. 

Of these, condition (ii) shows that U must be linear in both 
V and W. It is also a scalar. Therefore it must be of the form 
F, V + P.'W, where P.F' are two vector functions of the position 
of a point {i.e, functions of r), but independent of V and W. 

Taking then U=P.V + F'.W, condition (i) becomes 

VMP-V) + V2(P'.W)=0. 

at all points of the liquid. Or, since V and W are constant 
vectors so far as the operation of y 2 jg concerned, 

V.V’‘P+WV*F'=0. [See page 49.] 
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Or, again, since V and W are perfectly arbitrary, we have 

V^F=0 and V^F' = 0 ... (A) 

at all points of the liquid. 

Again F.V + F'.W being written for U, our condition (ii) 
becomes n. V(F.V-f F'.W)=:«. (V + Wx;*); 

or since V and W are arbitrary. 

7^.V(F.V) = ?^.V, 

and «.V(F^W) = ?^.W x r 

at all points r on the surface of the solid. If for a moment 
we write ^(Sr) and i^(8r) for SF and SF' respectively, we have 
by §30, p. 40, V (F.V) = </»'( V) and vF'.W)=^'(W), and 
therefore the above relations may be written 

n,Y=u4'{Y) = Y.<i>(n), 

and «.W X r = w.^'(W) = W.^(m) ; 

or, «=<A(m) = w.vF *) 

^ ( ... (B) 

and r X V F', ) 

since V, W are arbitrary. 

Thirdly, condition (iii) becomes 

V(F.V) = 0 and v(F'.W) = 0 
or il/'(Y) = 0.a.\\dtl/'(W)~0 ... (C) 

at infinity for all arbitrary vectors V and W. 

Conditions (A), (B) and (C) will uniquely determine the 
vectors F, F' and so U being determined, we have 

SiTi=/UvU. da. 

=/U( -f- W X r).day by the surface condition (ii), 

U being written for F.V-f F'.W. 
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4. To calculate Tg, the kinetic energy of the solid : — 

We have, 2 T 2 =f{Y + W x r). (V + W x r) dfn, dm being an 
element of mass of the solid at the point r and the integration 
extending through the mass of the moving solid. 

That is, 2 T 2 =J[V.Vh- 2 V.Wx ;* + W. (rxWxr)] dm [see 

p. 6]. 

~m V.V -\-%7n r. V X Wh- W.<i)(W), 

where m denotes the mass of the solid, r the vector to its 

centre of gravity (so that fr dm — m r), and w (W) has been 

written for Jr X (W X Writing this integral in the form 

/(r.rW— Wr.r)flf »2 [p. 6], we see that a> (W) is a self conjugate 

linear vector function of W. Clearly also <t)(W) represents 

what the angular momentum of the solid about 0 would have 

been, if O had been fixed; and W.o)(W) = constant, is, for 

variable W, the equation of the momental ellipsoid of the 
solid at O. 

5. We can now write down the expression for the kinetic 
energy of the system in terms of V and W, Thus, 

2T=2T,+2T2 

= JU V da + JU W X Tda + mYN + x W + W.a>(W). 

For R and G, then, we calculate V.T and V«T from this 
expression for T. Noting that and V 

if ^ is self conjugate [p. 28], and that V.U = F and V.U = F', 
we write down almost immediately 

2 R=2V,T=jUrfa + jPV.f?a+/PW xr.da + lniN + imW xr 

and 

2 G = 2v«T= jF'V.rffl-hJ \Jrxda-{-J F'W x r.^a + 2;;2r x V 

+ 2a)(W), 

6. There are now eerlain relations among the integrals 
occurring in these expressions for R and G. Just to obtain 
these we prove generally that if V is any constant vector, 

S S 

/PV. (^ffl.VP')=JP'.V(rfaVP), 
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E, F' being any two vector functions, satisfying v*F = 0 and 
V’F =0 at all points within the closed surface S over which the 
integrations are performed. 


If we put, as before, S F=(^(Sr) and 8F'=i/-(&'), then (fa.vF' 
and da.'^F=<t)((la). If then Cis any arbitrary constant 

vector 

S S 

C./FV. (r/«.vF')=/P.CV.,/r(^/fl) 


S 

=;F.C./r'(V).^^a 

S 

=/P.CV(F'V).^a. 


T 

which by Green’s Theorem=/v.[P.CV (F'-V)] (h, the volume 
integral being taken through the volume T enclosed by S. 
Using now (i) § 30, p. 38, and the relation v®F" = 0> 


S T 

C./FV. (fZa.vF')=/V(F.C). V(P'-V) dr-, 


S T 

Similarly, C./P'.V (^«. vF)=/V(F.C).v(P'.V)^/t. 
C being arbitrary, this proves our theorem. 


* 

In the application of this theorem to our problem, we have 
to note that the region of integration would be that between 
the surface of the solid and a sphere of infinite radius, and the 
question of convergence of the integrals would arise. This 
question has been discussed by Leathern in “Volume and Surface 
Integrals used in Physics,’* Sections IX and XI. 

7. Using this theorem now and remembering our surface 
condition (B), viz,, <f>(da)=da and \l/(da)=r x da, we have, among 
the integrals in the expressions for R and G in § 5 

JFV. ^a=jFV. <fi(da)=zf¥,Y<l>(da)=fF,Y da, ... (o) 


JFW xr. da = JFW. rxda z= JFW. ij/(da) = /F'.W <fi (da) 
^fW.Wda, ... (i5) 


• t • 
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i¥'Y.da=jF'Y 4 {da)=SF.Y^{,la)=!F.V,X(la, ... (y) 

s 

iWWxr.da = fVW.rxda = /F'W. (da)= P'AV ,p{da) 

=fF'.YVrxda ... (8) 

Hence those same expressions for R and G may be written 
R= jUda + mY + mYf Xr= [/P.Vr^a + w/V] + [/F'.Wf?a + mW x ?], 
G =/Ur xda + m'r x V + q)(W) = [/ F. V rxda-F ml- x V] 

+ [;F'.Wrx^i/ + <»(W)]. 

We write now R=<^iV + <^jW, where <^2 are the two 
linear vector functions, 

and <^,W=/P'.W^a + mWxr. 

Since the conjugate of «. r A is « r. b, (a) shews that </., is 
self conjugate. The conjugate of W 

=/F' W. da-mWxr 

= jP.Wr xda-\- mr X W, by (y), 

so that ’i>\Y=lY.Yrxda + mrxY. Hence we may write 

G=,^',V + 0,W, 

where <^,W=/P'.Wr x rf«+<o(W). Since <.(W) is self conjugate 

(p. 78 ), (8) shows that ^3 is self conjugate. 

Thus we may write 

R=<^iV + ^,W 

and Q=^',V + .#.,W. 
where are self conjugate functions. 

This fact alone,— of d,] being self conjugate and d>, and 

<#> , being conjugate— could of course be deduced directly from 
( 2 ), p. 74 .. 

8. Considering now the case where no extraneous forces are 
present, we have, putting C, X=0 in equations (1) of page 74. 
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cU 


= Rx W, 


and 


dt 


= RxV+GxW. 


The three well known integrals follow immediately. 

(i) Multiplying the first equation scalarly by R, we have 

R. =(RRW) = 0 

(it 


R.R = constant. 


(ii) G. +R.-^ =(GRW) + (RRV) + (RGW)=0 

at dt 

R. G= constant. 

That is, the pitch of the wrench (R,G) which is^^ is constant 

iv.il 

Again, if r = j which we know is the perp. from O 

R.R 

on the axis of the wrench (R,G), 


(!r=J_ riE. xG+Rx-^] 

(It R.R '-dt dt ^ 


= 1 [(RxW)xG + Rx(Rx V) + Rx(Gx W)], 

R.R 

which by the last formula of p. 6, reduces to 


[(RxG)xW + Rx(RxV)], 


R.R 


that is, -f- =rx W — Vj 

dt 

where V, = — ^ ^ = component of V perpendicular 

R.R 
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If then 



denote the rate of chanae of 





-rxWH-Y = 





where 



denotes the com 


ponent of Y parallel to R, i.e.^ to the axis of the wrench (R,G). 
It follows that this axis is a fixed line in space, its direction 
being obviously constant from the first of the equations of 
motion. 

But there was no special point deducing these results from 
the equations of motion, as the fact they express, viz.y the 
constancy of wrench (R,G) in ease no extraneous forces act, 
was obvious d prion from the theorem that “ the * impulse * of 
the motion (in Lord Kelvin’s sense) at time i differs from tlie 
impulse at time t ^ by the time-integral of the extraneous forces 
acting on the solid during the interval / — — [Lamb § 119] — 
of which theorem it is only an analytical ex])ression that we have 
in the equations of motion. 

(iii) V ■2T = R.Yh-G.W, 



Using (:l) p. 81, R.l^ +G. 


dN , ^ dW , . dR . da 

— -pG. reduces to V.-,- H-V>. — 
dt di dt dt 


=V. ^ +AV/^ =(YRW) + (WRV) 

di di ' dt ' ^ ^ 

+ (WGW) = 0 
/.Tss constant. 


9, These three integrals however are not sufficient to 
determine the motion completely. We require three more scalar 

9 
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(or one vector) integrals for determining the two vectors V and 
W. The difficulty of finding them is avoided in two cases : — 

(i) When R=:0, in which case we can utilise the known 
solution of Poinsot in Rigid Dynamics. 


(ii) Steady motion, when =0 and — =0, so that t 

(It dt 

is got rid of altogether from our equations of motion. 

10. The ease R=0, is fully worked out in Lamb*s Hydro- 
dynamics, § lid. It may be interesting, just by the way, to 
put the solution in vector form. 

The equations of motion reduce, when R = 0, to 


dt 


= Gx W. 


Again, putting R = 0 in equations (3) p. 81, 
we have 

Now (^3 is self conjugate. Also, the conjugate of a /product’ 
of linear vector functions being the ^ product^ of their conjugates 
taken in the opposite order (Gibbs’ Vector Analysis, Chap. V, 
p. 295), and further and therefore also being self 
conjugate, the conjugate of is itself. That is, G is a 

self conjugate linear vector function of W. 

Hence, G=-|. y.. (G.W) [see p. 28.] 

9 

It follows that Poinsot^s solution is applicable. 

II. But the case of steady motion does not seem to have 
received yet the attention it deserves. Two simple particular 
eases are well known — the permanent translation and the 
permanent screws that we have when R = 0. It is proposed to 
have a general investigation of the question here. 

If and ^.^-are both zero, =0 and =i0, andour 

dt dt dt 
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equations of motion of page 8*2 reduce to 

RxW = 0 

and R X V + G X = 0 

The first equation shows that R is pai-allel 

R=— .fW, where a? is a scalar. Substituting in 
we have 

(^'V-fG)xW = 0. 


... (I) 

to W, or 
the second, 


Again, putting R=— .rW in equations (3), p. 81, we have 

... (5) 

Therefore, - (a-V + G) = x* 61 . W + a-e, W + W = 12 W , 
where, (?, =<^7*, =<^7*.^, 

and 12 has been written for the linear vector function 


+.i^a 

We have, then, nWxW = 0; (11) 

and our conclusion is that in any case of steady motion W should 
he parallel to OW. 

12. By the rule for the conjugate of a product of linear 
vector functions quoted in § 10, we see immediately that 6,, 0,, 0, 
and therefore 12 also are self conjugate linear vector function! 

We recall now Hamilton’s theorem of the latent cubic of a 
linear vector functions— In general, for any linear vector 
function there are three vectors, say X,, X„ X,, of which the 
directions are left unaltered by the operation of that function. 
fftf ^3 are the roots of the cubic equation 

+ m'g—m = 0 , 

where ^[a-<^^X<^y 1 

a-/?xy apxy ' 


and m" = 

a,pxy 
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ay py y buing any three arbitrary noii-coplanar vectors, then, 

i>^ 2 =!/ 2^2 <^A 3 =:^ 3 X 3 , The cubic in y is 

Called the latent cubic of <l>, and its roots the latent roots, and 
the vectors A.^, A^ which retain their old directions after the 
operation of are called the axes of If </> is self conjugate, 
the latent roots are all real, and the three axes real and mutually 
perpendicular. [Killand and Tait’s Quaternions, Chaj^ X.]. 

Applying this theorem, we conclude that for any .c there are 
three mutually perpendicular directions for W, corresponding to 
any one of which Ave mav have a case of steady motion. 
Further, if y is tlie latent root of corresi^onding to the 
axis AV, f2Vy=yAV, and we have aV + G=-yW, or 

G=-.tV-^V. 


Thus our impulse is given by 

R=~.tAV 

and G=-.tA"-y\V. 



(7) 


L‘3. To construct a screw, therefore, such that the corres- 
ponding motion of the solid may be steady, we find an axis 
of the linear vector function n(r) corresponding to any si and 
take W along this axis. V then is given by (5), p. 85 to 


be — .r<^i ^AV — W. 


Draw the vec tor ^ an d through 

AV* W 


its extremity draw a line parallel to AV. This line is the 
axis of the screw. If then the motion is started by the 
impulse R= —.(AV and G= — 'A^—yAA'', the solid will continue 
to have the steady twist on our screw, the angular velocity 


being AV and the pitch 


AV.V 
AV • AV 


It is easily seen the axis of an}^ screw and the axis of the 
corresponding wrench coincide. For these axis being parallel 
to AV and R respectively are themselves parallel. Further, 

RxG _ -.:Wx(~.V-y/AA^)_A\^xV 
R.K (-a:AV).(-..W) AV.AV ' 


since, 


these lines are 
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drawn tlirough the same point, so that the lines coincide and 
either line is fixed in space. 

Pi P denotes the pitches ot' the screw and wrench 
(ll.G) respectively, 

R.R lA\).{-x^X) 

14'. We conclude then that any being specified, three 
mutually perpendicular (though not necessarily intersectino-^ 
screws of steady motion are determinate, except only as to the 
magnitude of W ; and again that, by varying ,r, the directions 
of screws of all possible steady motions would be obtained by 
solving the vector equation n(r) x ? =0, or 

xr=0, in other words, by finding an axis of the linear vector 
function ti. 

It would be interesting now to enquire to what values of x 
correspond the two cases mentioned in §11, t;>., permanent 
translation and steady motion with R = 0. 

Since, in any case, R= — -i W and W is not supposed to be 
infinite, R=0 would make .t=0. O now reduces to ^3, 

G = - W = - 6(3 W, V = - <#,7 * W . 

Since, again, when the motion of the solid is one of 

translation only, W = 0 and the screw reduces to V only; and 

since, in any case, R is supposed to be finite, we see that in this 
ease a; tends to infinity in such a manner that .i;W=— R. 
The equations (6) of p. 86, being written 

(..e,W + 0,W + i-6i,W) X R=0, 

.T 

we see that in this case 0,Rx R=0, i.e. <^7'Rx R=0. 

Hence the axis of the wrench, and therefore the axis of the 
screw also is parallel to r, whei'e r satisfies <l>~'rxr=0. That is, 
Y is paiallel to an axis of <^>7^ or to an axis of <^i, since 0, and 

axe co-axial. 
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These permanent translations, of course, could be worked 
out more directly by putting W = 0 in (4; of p. 85 and (3) of 

p. 81. 

15. 1£ we eliminate .t from 


we shall obtain the whole assemblage of lines, to one of which 

the axis of the screw corresponding to any ease of steady motion 
must be parallel. 

Multiplying sealarly by and respectively, we have 

^ I r .0, r X r + 0 , r.e, r X r = 0, 

and x?-H-e,r.e,rx/=0. 

■ _ e,r.(9,rxr 

' Lfl^r.0^rxrJ O^r.e^rxr 

or, {e^rfi^Txr) {O^r.O^r x r)-{e^r,e^r x ry , 

which is homogeneous and of the sixth degre*^ in the tensor 
of and represents therefore a cone of the sixth order, to some 
generator of which the axis of every steady screw must be 
parallel. 

Obviously, the axes of (he linear vector functions 6^ and ^s» 
— which we came across as giving the directions of screws for 
the special cases, § 14, — all lie on this cone, for the equation is 
identically satisfied if we ptit either $^rxr=^ or 0^ xr=0. 


16. The motion of the solid beiiiff thus a twist about a 

o 

screw of which the axis is fixed in space and pitch constant, 
it is almost obvious that each individual point of the solid would 
be describing a helix about the fixed axis of the screw and 
having the same pitch p (except of course points on the axis 
which would move along the axis). For, if referred to a fixed 


origin on the axis, the position of any point P of the solid 


is specified by r,(OT=;-), the velocity of P is^=/?W + Wxr. 

at 
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(Ir 

Jt 


) =WXA 


[Wx(Wxr)]=0 

which shows that (W xr).(W xr) is constant, or the magnitude 

of Wxr is constant. That is, the distance PN of P from the 
axis of W is constant. 


The formula for^^shows, moreover, that the velocity parallel 

to W is and velocity in plane perpendicular to W is W x r, 

that the motion in this plane is instantaneonsl 3 ^ a circle 
(of which the centre is N and radius NP) with angular 
e ocity | W | , Hence P moves in a helix of which the 

pitch = I _ „ 

I W I 

^ 17. For the maintenance of a motion of this type, appro- 
priate forces must be continuously acting on the solid, for we 
now that the only motion a solid is capable of under no forces 
18 either one of uniform translation, or a uniform translation 
Mmbined with a motion of rotation about a principal axis of 
e solid at the centre of gravity. We shall just verify that 
111 presssures exert on the solid just the force and couple 

necessaiy for the maintenance of the sort of motion that we 
have here. 

Considering the general ease (where the motion is not 
n^ssarily steady), let f', denote the force and couple which 
t e fluid pressures on the solid are equivalent to. The linrar 
and angular moments of the solid are respectivelj’ 

B,,=«(V+WxT‘) and G,=»ir x V-(-<o(W), 


• • » 


( 8 ) 
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these being just the terms of R and G that are obtained from 
by the oi^eration of and V». respectively. For the solid 
alone, therefore, the equations of motion are 


(IK. 

i 

~Tf 

r/G 


= R, xW-ff 


and =R, xV-i-G, xW + A'. 


If we farther put now R, and Gi ==/Ur x we have 

R = Rj-fRjj, G = G^+G 2 . From our equations of motion 
therefore of p. 32, we have 


xW-f aucl 


(It 


rIG 

(If 


= R, xV + G,xW— V 


Hence, +R, xM^ and 

(It (If 

+ R^ X V + G, xW. 

dK ^ 

Considerinar steadv motion now, for which — =0 and 

” ' dt 

^^ = 0, we have 
dt 

$' = K^xW and X'=:R, x V + G, xW, 
or, using (4-) p. 85. 

f =— R^ X AV and -R, x V-G, x W. 

These formulae, for the special case of steady motion, could 

of course be obtained directly by piitting^^- =0 and ^ ^ 


in the equations of motion of the solid above. 

Substituting now the values of Rg, Ga given in 

* 

obtain after slight reductions, 

^'=W x ?//(V 4- W Xr) 
and A^ = w/(VxW)xr — (u(W)x^'. 


(8). 


we 
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If, again, we refer the motion to the centre of gravity of 
the body as origin, so that r=0, we have 


I' = W X V 


and X' = Wxo>(W), 

where (t>(W) now is the an 
its centre of gravity. 


(9) 


gular momentum of the solid about 


We notice that and both vanish only if either 
(i) W = 0, or (ii) V, W and a>(W) are collinear. The first is 
the ease of one of the three permanent translations. In the 
second ease, since W and w(W) have the same direction, W is 
along an axis of the linear vector function w, i.e.j along a 
piineipal axis of the solid at its centre of gravity. The axis 
of the screw therefore is a principal axis of the solid at its 
centie of gravity. These two, of course, are precisely the cases 
which we expected d prioii that f, X' should vanish. 


We show now that iu the general case of formula 

(9) are just the force and couple that would make >he solid 
continue to have its screw motion represented by V, W at the 
centre of gravity. Since the motion of the centre of gravity and 
the rotation of the solid about the axis W at its centre of gravity 
can be considered independently, we show that the velocity V 
of theC.G. is maintained by and then, regarding now the C.G. 
at rest, that the rotation W by itself is maintained by X' ; or, 
what comes to the same thing, that the mass-acceleration of the 
C.G. 18 equal to and that the rate of change of angular 
tnomentum about c.5. is equal to X'. 


For, if referred to a fixed origin O' on the fixed axis of the 
screw, the position of the C.G. is specified by r, then its velocity 



=sj9W 4*W X r, 
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7n^ X V, which is our . 

Again, the angular momentum of the solid about the C.G. 
is w(W), and in calculating its rate of change we take the C.G., 
which is now' supposed to be at rest, as our origin. Since the 
'^origin system” [see Silberstein’s foot note cited ' on p. 74] is 
fixed in the body and rotates with it with angular velocity W, 

the rate of change in question = — aj(W) x W 

= Wxo.(W), [V =0], 

which is our X'. 

18. Having thus considered the general character of steady 
motion in the preceding articles 11-17, we would next turn our 
attention to the question of stability of these steady motions. 
Before, considering this question, however, it would be convenient 
to summarise here a few propositions on the linear vector function 
which we shall presently have occasion to use. The more 
important ones are taken directly from Kellaud and Tait^s 
Quaternions, Chap. X. 

(i) Prom the well-known relation 

r{aPy) = (r^y)a -f (rya)/3 + 

it follows easily enough, — say, by writing in the form 
Xr-\-a.rb-^cx7’ of page 24 — that 
{aPy)4>{r)-{Tpy)tf>{a) + 

Here a.,P,y are any three arbitrary non-coplanar vectors, and r 
any fourth vector, and ff> denotes a linear vector function, (a^y) 
etc. of course, as explained on p. 6 stand for a.fi x y, etc. 

(ii) For inverting the functign ^(r) we have, if we denote 
the inverse function by 

» (X x /i) = «^'X X 0 
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where = and A, ft any two vectors. 


If we introduce the facction ^ by the definition, 

we may write symbolically = 1 //. 

(iii) lor inverting or where ^ is a constant 

scalar, we have 

{m + rn'g + ■’)(<#, +y) - (X x m) = +^)X x (</,' + ^)/x 

= i'l'+ffX+ff'‘) (^ xm ), 

where m, m' , m" have the values defined on p, 86, viz. 

m{aPy) = = 2(a<^^<^y), >«"(a/Sy) = 2(;8y<^a), and 

the function x is defined by 

x{^^t^) = <l>'kx fj.-\-Xx<l>'fx, ^ 

(iv) Integrating over the surface of a parallelopiped of 
which the edges are the vectors o,/J,y we get easily enough 
i<l>r.d<T^'S,{py^a). 


Hence our 

(w) 


=scalar constant of 



(v) We have also for any vector X. 

j/i"X=:<^X -f- ;^Xj 

or, symbollieally, W 2 "=<^+X. 

(vi) We already enunciated Hamilton’s theorem of the 
latent cubic on page 86. We only note here that if for any 
function jw=: 0, the product of the three roots and y, 
of the latent cubic vanishes and therefore one of the roots, say 
y,, is zero. It follows that for such a function there exists a 
vector Xj. such that <#.X,=0. Conversely, if for a function 4> we 

can find a vector Xj such that w for that function must 

vanish 
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Keeping* to Gibbs’ notation, we shall speak of vi for any 
function 0 as its (Jclerminant and denote it by | <^ | . 

(vii) About the determinants of linear vector functions, 
we have the theorems that the determinant of a ^product’ of 
linear vector functions is the product of their determinants, and 
that the determinant of the ‘ quotient ’ of two linear vector 
functions is the quotient of their determinants. Thus the 

determinant of I I = 

I I 

are the determinants of respectively. [See Gibbs, 

p. 

(viii) The inverse of the ^ product ’ of any number of 
linear vector functions is the ^ product ’ of their inverses taken 
in the opposite order. [Gibbs, p. 293]. Thus the inverse 

19. To examine now the stability of any particular mode 
of steady motion, we consider as usual the effect of a small 
distuibance given to the sj^stem. V, W being the linear and 
angular velocities of the solid for the steady motion, and Cr 
the corresponding impulse, if SV, SW are the additional linear 
and angular velocities imparted to the solid by the disturbance, 
the total impulse that .vould generate from rest the new 
disturbed motion is represented hy 

R + 8 R=<^,(V + 8 V) + ^,(W-f 8 W) 

and G + 8 G=<^\(V + SV) + <^ 3 (W-fSW), 
so that 8 R=<^, ( 8 V) + 1^3 ( 8 W) 

and 8 G=<^' 3 ( 8 V)+</) 3 ( 8 W). 

Our equations of motion of p. 82 now are 

(R + 8 E)= (R + 8 R) X W + SW) 

fti 

and^ (G + SG) = (R + SR) X (V + 8 V) + (G + 8 G) X (W + 8 W), 

at 
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which, by the condition of steady motion [§ 11, p. 85], 
reduce to 


dZR 

dt 


= Rx8W + 8Rx'VV, 


and i5^=(RxSV + 8RxV + Gx8W+8GxW, 
dt 

if we neglect the terms 8Rx8'W, 8Rx8V and 8G x 8W. 
Using (7) of p. 87, these equations further reduce to 


d8R 


= (8R + 8W) X W, 


<f8G 


( 10 ) 


and^ =(8R + a8W) X V + (8G + .. 8V + y8W)^W. I 

dt j 

20. Let us put now SV=t'*e’*' and 8W=M;e’'', w^here Vytv 
are two vectors independent of ty so that 

8R=RV' and 8G=GV', 

where R'=<#»jV+</»,w ■) 

[ ... ... (11) 

and ; ' 

and ^ =ne"'R', ^ =ne"'G’. 
at at 


Equations (10) then are identically satisfied provided 


or, 


nR' = ( R' + .1 lo) X W, 

and nG*=(R^ + .i'w) X V + (G* + .it’ + yt(;) x W 

+ -N 

V * 

and ) 


where ^iVX W, 

+aj)w X W, 

*»v=n4,',v->f>^vxY-(^',+x)vxW, 
and **w=n.^,u>-(<^, +*)«> x V-(t^, +y)w x W. 


I ... (12) 


... (13) 


» • • 


( 14 ) 
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irom (18) noTV we have r=— 

and ~$ 3 ^ 7’<^2 7r4-4>j?r=0 


(15) 

(16) 


It follows that the latent cubic of the function— 4 > 3 <E> 7 * 4>2 
has a zero I’oot and w is the corresponding axis. Thus the 
determinant of the function vanishes and we have 

[See (ii) and (vi), § 18] 

“j fit y being any set of three non-coplanar vectors ; that is, 

-2[4>^a^>^^X^>34>7>4>,y] + 2[d>*a.$3$7>4.,/?x$3$7'^>,y]=0. 
Dividing out hj' (afiy), and noting that by § 18, (vii) p. 99, the 


determinant of $ 3 $ 7 >$ 3 = — I I ^2 ., ive have 

*1 I 




(“^y) 


2[<E>4a.5>,3^x4>3$7’$ay] 


i * * 

i t i 




This relation (17) gives us an equation for the appropriate 
values of n, and only when the roofs are all imaginary, would 
the motion be stable. 

21. Summary of results to be proved . — It is very tedious 
n'ow to calculate all the co-efficients of the several powers of n 
in this equation, and even when these co-efficients are calculated, 
it is found impossible to come to any definite conclusion as to 
the conditions of stability. The solution of the problem with 
any degree of completeness was therefore given up as hopeless 
by this method, and there is no point going here over all 'the 
elaborate analysis for obtaining the full equation in n. What 
we propose to do is just to obtain a few special results which 
can be proved without much trouble. Thus we shall ^how, m 
the first instance, that in the most general case our equation 
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for H is of the seventh degree, and that the co-efficient of «' is 
determined solely by the nature of the solid and does not vanish, 
unless some special limitation is imposed on its shape, say, by 
way of symmetry. If will follow that, being of an odd 
degree, this equation must have at least one real root, and that 
in general therefore, for any arbitrary ,r, we shall not have any 
stable steady motion at all. We shall next show that for some 
special values of r, the term independent of a in our eipiation 
vanishes, and then neglecting the root ?i = 0,- of which the 
effect is simply to add a constant vector each to the values of 
8V and 8W, and which therefore does not affect the question of 
stability either way,— we shall have our equation reduced to 
one of the sixth degree ; and since an equation of an even 
degree may have all its roots imaginary, it is possible now for 
the motion to be stable. Instead of groping about then foi the 
stable screws among the infinite system of steady screws, we 
shall be sure of this one fact that if any stable screws are there, 
they would be found only among the group determined by 
those values of x which make the term independent of a 
vanish. It will be shown again that in general there are six 
such values of x, and our conclusion would be that in any case 
there cannot be moie than 18 (=6x3) stable steady screws. 
We shall show finally that .<;=0 is one of these six values of x, 
and we shall conclude by working out this ease completely. 

22. To obtain the degree in n of the left hand side of (17), 
let us first express R' iu terms of w from the first of the equa- 
tions (12) of page 101. 

Writing ,^R' for «R'-R' x W for the time being, this equation 

IS <;>R'=nR'— R'x W=a:icx W. Hence by (ii) of §18. p. 97 
we have ’ ^ ’ 

W) =x<l>hv x <^'W, 

where f mis the conjugate of 4,w and is therefore rnv + wxW 
and and 


M«^r) = («a-a X W).(nfi-pxW) X („y-y x W), which. 

wf r to be equal to 

n®(a^y) + nW.W(a^y), bo that w=M(n* + W. W). 
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Thus we have, + W.W)R'=,i-(MM;+tox W) xmW, 
or, (m» 4 . W, W)R'=a;[«(!tj x W) + (w x W) x W] 

=,.[m(wx W)+W.j(,-W— W.W m!], p. 6 , ... (18) 

23. Now from (11) p. 101, B,'=<|>^v+(j>^tv. Therefore, 

’R'. Equating this value of v and the value 
in (15), p. 101, we have, hv using (18), 


— r=4>, 


+ w W X W) + W. w ^7 > W 




~ n»4W.W 

+ W.W<^7 ' (<^ a + .r)lO— W.w.t<^7 • WJ 
Looking up the form of ^3 now in (14), p. 101, we see that we 
can write w= L— — fn® A + n*B+nC + D]tt;, where 

w®+W.W 

A, B, C, D denote certain linear vector functions which do not 
involve n in their constitution. We have, for instance. 


and again, it is easily seen also, that 
— Dw=^i[W.W<^7'(<^2 +.»■)«;— W.m?.i^7^WJ x V 

+ (<^'2 + .r)[W.W^7^(<^3 + 

=W.W[(i ^2 + xY + (<t>\ + u;) 07‘(<^8 + .p)w X W] 

— W.^ — V + ( <^' j + a :) ^7 ^ W] X W 
= W.W [(<^3 4- .c)w X V + (^3 4- n)w; x W] 

-a!W.«;[2u^, W 4 -^,W] x W, by (5), p. 85, 

= W.W[(<^g +x)w X V4- (< 3 (»a +n ^ 4 * W.-wF) x WJ, ... (20) 

where fl, 0^, 6^ are the lineai* vector functions defined on 
p. 86 and F has been written for the vector function of W, 

-^[2^0.W+W 
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24. Our of (17), p. 102, then, is of 

the form 

where nA' + B' has been written for so that 

Aa=;<;^ja, and — B'o. = (<j>^ + .i;)a xV+(<^ 3 + 2 /)ax W. 

Obviously A', also do not involve n. We may write therefore 


2[0*a.‘3>*/3x$3<l>7^4»,y] 


i^Py) 


n* + W.W 


[a5n*+a*n‘+a3n3+a,^n*+ain + ao] 


where the o s are all independent of n. The values of and Uq 
niay be written down immediately. Thus 

(a/Jy)a.=2[A'a.A'^xAy]=2[<^,a.^,;8x<#>',</.7'<^,-),], ... (21) 


and 


(a/3y)a„ = 2[B'a.B'/JxDy] 
25. Again, with the same notation, 


... ( 22 ) 




_ 1 1 

(^) ( n*+W.W) » 2t("^' + B')“. (n»A + n*B + nO + D)/J 
' # 

X (n*A + n*B + nC + D)y], 

which is of form K,n' +»',»" +a +a'„], 

where also the co-efficients of the powers of n are independent of 

n. Clearly also, 

(a^y)a',=2[A'a. A)3 X Ay]=2[<#»*a. 


... (23) 

(a/3y)a\=5[B'a. D^xDy] (24) 


and 
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26. We try to form an idea now of the other terms 

I I » I ^2 I I ^3 I that occur in (17), p, 102, 

(i) We note in the first instance that </>(r) being any liner 
vector function and a any constant vector, the determinant of 
<f>r X a vanishes. For, 

(«^aXa). (<I>P X a) X ((f>y X a) =[<t>a X a]. a(^y3^'ya)=0, [page 6.] 

It follows that the parts independent of « in and viz.^ 

— 0jrxW and — (^j+a!)rxW have their determinants equal to 
zero. 

(it) We note again that if are any two linear vector 

functions, the determinant of is easily found by writing 

out the expanded form of 

(^) Wiy+'^«y)l> 


to be where Wi,Wg are the determinants of 

respectively and /x, /t' are determined from 

/x(a^y)=5[^ia. 

and M'(«^y)=2[<^ia. <^a^x<^,y] 

It follows that the terms independent of n vanish in both 
I I and I I . We find, in fact, by working out 

[(n<^iO~<^jaX W). .^,/?xW)x(n<^,y— <^iyxW)], 


as in § 22, that 


( =wijw(n*+W.W), 


where is the determinant of 
We may write also 


I I =Wan® fftjW^+Cati, 

I *3 I =Ws«®+65n»+C3« + da, 
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where is the determinant of either of the two conjugate 
functions </> 2 , <l>\ and m 3 is the determinant of <^ 3 , and the b't, cs 
and <rs are all independent of n and in general none of them 
Tanishes. 

27. For our equation in n then, corresponding to (17) p, 102, 

we have 


mi(n* +W.W) 

+ + 

n*+W.W 


+ 4- ». + a' o 

(n® +W.W)* 


••• ... ••• ••• (25) 

Multiplying out by (w*+W.W)®, we get an equation of the 
seventh degree in n. 

The co-efficient of n’ 




1 

i^Py) 




' I I , 

by (21) p. 106 and (23) p. 107, and (vii), p. 99. 

Remembering now theorem (ii), § 26, we see that this co-efficient 

= I ! =— I ^3 I ,^3 having the same definition 

as on p. 86 . Obviously, our co-efficient of n’ does not vanish, 
unless the functions and <t>^ are given in special ways 

unless, that is to say, some special restriction is imposed on 
the shape of the moving solid, for these functions are solely 
determined by the form of the bounding surface of the solid. 
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28. We may wiite down also the term independent of n in 
the same equation. This term directly 

= (W.W)^rf,-a,W.W+a',-W.W 

Now if, as in § 24, we write n A'+B' for by (ii) § 26 is 

the determinant of B'. That is, 

d,(afiy)=i(B'aB'/3B'y), 

« 

The values of a„, a' „ are given by (22) and (24), viz., 

(a/3y) a.=2[B'a. B';8 X Dy] =W.W2[B'a.B'/? X D'y], 

(a/3y)a'„=2[B'a. DyS X Dy] = ( W. W)«2[B'a. D'/3xD'y], 

if a new function D' is introduced by tlie definition 

-DV= Dr, which by (20) p. 106 


= + X V + [(^3 W.rP] X W, 

where F=— [2,r^, W+6',W] 



Hence, (W.W)®d^— a«,W.W+a^ 


= [(B'aB'^B'y-2(B'aB'j8D y) + 2(B'oD'ySD'y)] 

=(W.W)“[ I B'-D' I + I U' I ], by (ii), § 26, 


But (B'-D>=[(n-y)r + W.rF] X W, 
and, therefore, by (i) § 26, | B'— D' | =0. 

Therefore, (W.W)*d^— a„W.W+a\=(W,W)* | D' | . 

4 

Hence our term independent of w 


=W.W[W.W I D' I 
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^9. tor further simplification of tiiis tei'm, we calculate 
generally the determinant of 

<^iyxv+<^,yxw. 

'#‘ 1 ) <#>a being any two linear vector functions. 

Since the determinants of both x V and x W are zero 

by (i) § ^6, the determinant in question is by (ii) § 26. 

~ (^) V)x(.^,yx W) 

+ 2(<^,axV). (<#>jj8xW)x(<^,yx\V)] 

Now 2(</.,axV)x(<#,,/3xV).(,^,yxW) 

=2(<^,a<^,/3V)V. (.^,yx\V) 

= WxV. 2[(<#..a./,,/JV).^,y] 

= WxV. by (i), § Ig^ p 9(5 

=WxV. <^,<^7'V' »i,(a^y), TO, being | <^, | . 

Again, 2(<^.a X V). X W) X (,^,y X W) 

= 2(<^,axV).W(iji,a<^,yW) 

=VxW. 2[(^,/S<;..yW)0,a] 

= VxW, '#>,<^I'W(<^,0<^,(8^,y) 

=VxW. m,(a^y), TO, being ' <^, | 

The determinant of <#>,»• xV+<^,rxW is therefore 

=WxV. [to,<^.,^7IV-to,.^,</,7‘W] ... 

=WxV. <^,i/',W],in the notation of (ii), p. 97 

18. 

30. Applying this result now, we easily calculate d, and 
I D I , 
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Thus, if n e write $ 3 , is by (ii), § 26 the deter- 

minant of B". That is, 

— d 3 =determinant of V + ((/)'a-f-.y)?*X W 

=WxV. [(0'a -h-O'Ai V-- 0 i(i/r' 2 -|-r*)'W'], by 

(hi) § 18, p. 97 

Again,— | D' | = determinant of (<^a +Jj)r X V'+/rx'W, 

[(26), p. 110] 

if/ris written for the linear vector function (^j +0)r-f W.rF, 
Hence, we have— | D' | =WxV. [/(^,+.c ir8+.c“)V 

where g denotes the ^ function for /, that is, 

g {aX^) ^faxf^ = [^3+fi) a + W. Fa] X [«^s ^ 

+ W. Fy9]. 

31. For the expression (27), only remains to be calcu- 
lated. This is done directly. Thus by {ii), §. 26, 

(a^y) = [(^a + i*!) ^x'W')]x[(^g+a!) yxW] 

= 2</»aa. W [(</>, +j;) /3x (<^a+‘») 7- W] 

= 2</»8a. W [Cc»-f U;x'a+<A\) ^Xy)]. W, [see (tit), p. 97] 

= W [.X® ()3yW) +. 1 ; (XsW^y) + (^jW^y)] 

= W. [.1!®2 (^yW) </»ga-f.T2 (iTaW^y)^, a 

-fS(V^,W^y) ^a®] 

which by (2), §18 = (a/ 8 y) W. W + .r<^2X2 W + «2 W], 

since by (ii), §18, </» 2 ^ 2 W = Wg W. 

/• C2=ifc® W. ^2 W -f .« W. ^ 2 X 2 W + W, W. 

32. The term independent of n, then, in our equation for 
n is, as given in (27), 

W. W [W. W I D' I - 1, 

7^1 J 
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where c^, and | D' | have the values caleulatcd in the last 
two articles. As explained in §21, it is necessary for stability 
that this term should vanish. A necessary, - of course, not 
sufficient - condition therefore that any one of the three steady 

screws corresponding to .c should be stable, is that a: should 


satisfy 


w. W I D' I - = 0 


» • 


... (28) 


Since now x occurs in the first degree in V [(5), p. 85], in 

the second degree in fi [p. 86] and P [(26), p. 110] and in’ the 

fourth degree in g, it appears from an inspection of the values 

of c„ d^ and I D' | given in the last two articles, that (28) 
represents an equation of the sixth degree in ,«■, and will, in 

gewml, therefore determine only six values of .r. No other 
than these six can possibly determine a stable steady screw. 

33. We prove now that ^=0 is always one of this set of 

SIX va ues of .r. It is only necessai-y to show that (28) is satis- 

neawhenaj = 0 . 

For, when a:=0, the expression for -DV in ( 20 ), p. no 

1C68 to 


reduces 




P ^ing zero, and 4,^+0 reducing to 4 > 3 + 03 , i.e. to 
m this case. Hence now by {Hi), p. 11 2 ^ 

- I D' I =WxV. [m. I 4>'AV’t> 


But, by (®m) and (tiw), §18, 




14, 




m 

I 4>',rr^4>, I = and (4>\4>V<l>,y^ =<l>V4>r4>'V. 


Hence, - I D' I = WxV. [m, 4 ,\ 4 ,-^Y^ 


m 




= .^WxV.[^<,^^V-.^.^,W],by(ff), 


§18. 
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Also, putting a; — 0 in the values of Cg, in §§30, 31, we 
have now W. W 

and -^ 3 =Wx V. [0', 


It follows that W. W I D' 


=0 identically 


when 


a: = 0. 

/ 

34. Thus the case a; = 0 satisfies our necessary condition of 
stability. It is a relevant enquiry then if any one of three 
screws corresponding to this case is stable. This is the ease 
(§14) when the impulse producing the original steady motion 
reduces to a couple alone, G= — = — yW ; R=0. 


For equations (12) of p. 101, we write now 

;2R'=R'xW 

ttG'=R'xV+ (G'+yzi;) xW. 


From the first, R' = 0 (for no vector can be perpendicular to 
itself), and therefore the second reduces to 

7iG'=(G'H-yzi?) X W ... 


Putting again R^ = 0 in (11) p. 101, we have 
V = — G' = — 

Writing for 0^ for convenience, we have G'==— 

and /, — — ( — ^ 4 *y) wx W = 0 

or — (<^— y)?/;x W=:0. ... (30) 

Putting the determinant of the function on the left hand 
side of (SO) to zero, we shall have our equation for the appro- 
priate values of n in this ease. 


The determinant in question is, by («) § 26 
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where, —d— determinant of ( — y) wx W = 0, by (/) ^ 26 ; 

///= determinant of <^ = | ^ | ; 

0;(W + „AV. W, just as in the 

calculation of t-g in § ;il ; 

and -'j{a^y) =^^ax<l>^). [,^_^)yx\V], which we lind to 

be zero when we write it out and remember also that <#> is self 
conjugate ; so that 6 = 0. 

Hence, our equation for n now is 
/« /r3+« [y^W. <^W-j,W. .)ixW + y«W. WJ=0 ... (til) 

[We may note, by the way, here that the co-effieient of 

m tins equation, ,dz., ,n or | , | , is tl.e same (but for the 

sign) as the co-efficient of in the equation for u in the o-eneral 

ease (see § -17, p. 109). Since this last co-efficient in any case 

IS independent of jr, and since equation (31) is only the reduced 

form, when '=0, of the general seventh degree equation, we 

see that when : =0, not only does the term independent of // 

vanish in the general equation, (as we have proved in § 33) 

but the co-efficients of all terms from downwards vanish 
too]. 

33. Now one root of equation (31) is n = 0, and the other 
two are j^iven bv 


= - L 
/// 


[y'-iw. 4,W-,/W. <l>xW + fnW. \V] 


since .t.W = //\V, and by (V) v' 18, this 

W. W 


simplilies to 


== - 


m 




The motion is stable, tl.erefore, if y is so chosen that 

1 


/// 




$ f • 




ispositiv. Now, by Hamilton’s theorem of the latent cubic 

U, p. 86) ,// may have any one of three values ’ 

say y 1 . y n , ^3 . of the equation ’ ’’ 
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so that /// =:^j -f = 
Putting, then, ij equal to 
expi’essiori (32) becomes 


!Ji!h. and t/i-- 

any one of these roots, 




1 

//i Ht.U 


i</i +I/i + >/s) ■ 3], 


I 2 ~ 

~rir^ (// 2 +//s) + 

.72 Hz 


IX. 


1 


y>i ^3 


(y 1 y^) (y 1 y.'j) ■ ■ • 


(33) 


The steady screw, therefore, parallel to that axis of which 
corresponds to tlie latent root 1 / ^ is stable, if the expression (33) 
is positive. 

Since, R being zero, the energy of tlie steady motion for 
any y is v U. AV = — ^ ij AV. AA", which must in any ease be 
positive, it follows that all the are negative. The expression 
(33) therefore is jiositive, if is numerically either the greatest 
or the least of the three numbers yj, ^ 3 , and it is negative, 
if is, intermediate in magnitude between 2 ya* 
last ease, therefore, the motion is unstable, and it is stable in 
either of the two other eases. 

AA^e may ])ul the conclusion in another form. Since — 
is always jiositive, — where k is a positive constant 
represents an eIIi[)Soid, of which the principal axes are in the 
directions of the axes of the linear vector function </», and the 
magnitudes of these princijtal axes are inversely proportional to 
y —y/j, s/^^hi \/^ 3 - Hence, the two steady motions 
for which the screws are parallel to the greatest and least axes 
of this ellipsoid are stable, and that steady motiop for which 
the screw is parallel to the mean axis is unstable. 
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